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PREFACE 


In  October  1981,  the  U.  S.  Army  Engineer  Waterways  Experiment  Station 

(WES)  was  requested  by  the  U.  S.  Army  Engineer  District,  Norfolk,  to  perform  a 

hydraulic  model  investigation  of  possible  hydrodynamic  changes  to  Chesapeake 
Bay  as  a  result  of  the  proposed  deepening  of  the  Norfolk  approach  channels. 

The  study  was  conducted  by  personnel  of  the  Hydraulics  Laboratory,  WES, 
and  its  subcontractor  Acres  American,  Inc.,  under  the  general  direction  of 
Messrs.  H.  B.  Simmons,  Chief  of  the  Hydraulics  Laboratory,  F.  A.  Herrmann,  Jr. 
Assistant  Chief  of  the  Hydraulics  Laboratory,  R.  A.  Sager,  Chief  of  the  Estu¬ 
aries  Division,  W.  M.  Dyok  of  Acres  American,  Inc.,  and  Dr.  R.  B.  Taylor  of 

Tetra  Tech,  Inc.  Testing  was  conducted  under  the  supervision  of  Messrs.  R.  0. 

Bruno,  Chief  of  the  Chesapeake  Bay  Model  Branch  (WES),  and  S.  R.  Rives  of 
Acres  American,  Inc.  Data  analysis  and  final  report  preparation  were  con¬ 
ducted  under  the  supervision  of  Messrs.  Bruno  of  WES  and  J.  R.  Pagenkopf  of 
Tetra  Tech,  Inc.  Project  Engineers  for  the  model  study  were  Messrs.  D.  R. 
Richards  for  WES  and  M.  R.  Morton  for  Acres  American/Tetra  Tech.  Additional 
key  personnel  for  WES  involved  in  the  model  study  included  Messrs.  A.  W.  Crunk 
M.  A.  Granat,  and  Ms.  V.  R.  Pankow.  Key  personnel  for  Acres  American  and 
Tetra  Tech  included  Messrs.  D.  G.  Dionne,  W.  E.  Hayes,  P.  S.  Jayne,  P.  A. 
Waltz,  and  Ms.  M.  F.  Capriotti.  Special  acknowledgment  is  made  to  Mr.  Crunk 
for  providing  computer  graphics  assistance  which  was  used  extensively  in  this 
report.  This  report  was  prepared  by  Messrs.  Richards  and  Morton. 

Commander  and  Director  of  WES  during  the  conduct  of  this  study  and  the 
preparation  and  publication  of  this  report  was  COL  Tilford  C.  Creel,  CE. 
Technical  Director  was  Mr.  F.  R.  Brown. 
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CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (SI) 
UNITS  OF  MEASUREMENT 


U.  S.  customary  units  of  measurement  used  in  this  report  can  be  converted  to 
metric  (SI)  units  as  follows: 


Multiply 


1 i 


To  Obtain 


acres 

cubic  feet  per  second 
feet 

feet  per  second 
inches 

miles  (U.  S.  statute) 
square  feet 


0.4G47 

0.02831685 

0.3048 

0.3048 

25.4 

1.609344 

0.09290304 


hectares 

cubic  metres  per  second 
metres 

metres  per  second 
millimetres 
kilometres 
square  metres 
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NORFOLK  HARBOR  AND  CHANNELS  DEEPENING  STUDY 


PHYSICAL  MODEL  RESULTS 

Chesapeake  Bay  Hydraulic  Model  Investigation 


PART  I:  INTRODUCTION 

Chesapeake  Bay 

1.  Chesapeake  Bay  with  its  tributary  estuaries  forms  the  largest  estu¬ 
arine  system  in  North  America.  The  190-mi le‘--long  estuary  varies  in  width 
from  4  to  30  miles  with  an  average  depth  of  28  ft  (Figure  1).  The  mean  annual 
discharge  of  its  126  freshwater  tributaries  is  approximately  70,000  cfs,  al¬ 
most  90  percent  of  which  is  contributed  by  the  Susquehanna,  Potomac,  Rappahan¬ 
nock,  York,  and  James  River  basins.  The  Atlantic  Ocean  provides  salt  water  tc 
the  bay,  producing  large  salinity  variations  within  its  boundaries.  The 
eastern  shore  is  generally  saltier  than  the  western  shore,  attributed  in  part 
to  the  predominance  of  freshwater  flow  from  the  western  shore  tributaries  and 
to  the  counterclockwise  tendency  of  flow  resulting  from  Coriolis  force. 

2.  Chesapeake  Bay  is  classified  geologically  as  a  drowned  river  valley 
estuary.  The  Holocene  sea-level  rise  inundated  the  Susquehanna  River  Valley 
to  form  the  bay.  Sedimentation  from  the  tributaries  as  well  as  erosion  of  the 
banks  has  contributed  to  maintaining  the  bay's  broad,  shallow  character.  The 
bay  is  classified  as  a  partially  mixed  estuary,  although  various  stages  of 
freshwater  discharge  and  tidal  and  wind  mixing  cause  portions  to  alternate  be¬ 
tween  well  mixed  and  highly  stratified.  Tides  are  semidiurnal  with  mean 
ranges  from  1  to  2  ft.  The  length  of  Chesapeake  Bay  is  such  that  a  complete 
tidal  wave  is  contained  within  its  limits  at  all  times.  Wind-generated  waves 
are  generally  less  than  3  ft  in  height,  but  larger  waves  can  occur  during  high 
wind  conditions.  Average  maximum  velocities  for  tide  and  wind-driven  currents 
range  from  0.5  to  3  fps. 

3.  Chesapeake  Bay  contains  a  variety  and  abundance  of  wildlife  that 
make  it  unique  as  a  biological  environment.  Its  estuarine  waters  are  rich 


*  A  table  of  factors  for  converting  U.  S.  customary  units  of  measurements  to 
metric  (SI)  units  is  presented  on  page  3. 
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with  nutrients  and  organic  material  that  support  some  of  the  largest  crops  of 
oysters  and  clams  in  the  world.  The  lower  salinity  portions  of  the  estuary 
provide  the  spawning  grounds  for  a  variety  of  fishes  including  the  striped 
bass,  shad,  and  herring.  Most  of  the  organisms  native  to  Chesapeake  Bay  are 
tolerant  of  natural  fluctuations  in  the  saltwater  distribution,  but  sometimes 
permanent  changes  provide  stresses  that  they  cannot  survive.  The  adverse  ef¬ 
fects  of  these  sometimes  man-induced  changes  are  often  realized  through  the 
depletion  of  existing  spawning  areas  for  fishes  or  by  providing  more  conducive 
environments  to  the  organisms'  predators  as  is  most  noticeable  in  shellfish. 
Since  man  is  capable  of  making  permanent  changes  to  the  bay’s  salinity  and 
hence  its  biota  through  his  navigational  improvement  efforts,  it  is  important 
to  know  the  possible  impacts  on  the  environment  prior  to  implementation  of  a 
project  so  that  reasonable  alternatives  can  be  explored  or  trade-offs  assessed 
in  a  responsible  manner. 


Norfolk  Harbor 


4.  Navigational  uses  of  Chesapeake  Bay  in  the  Norfolk  area  are  of  great 
importance  to  the  Nation  and  the  local  communities.  Due  to  its  naturally  pro¬ 
tected  harbors,  the  Norfolk  area  has  historically  been  the  home  port  of  naval 
activities  since  colonial  times.  Commercially,  Norfolk  has  played  a  major  role 
in  east  coast  bulk  shipping  for  many  years.  Its  closeness  to  the  Appalachian 
coal  fields  and  connecting  rail  lines  has  helped  it  become  the  largest  coal 
exporting  port  in  the  United  States.  However,  with  the  current  trends  toward 
deeper  draft  bulk  cargo  vessels  and  an  ever-increasing  demand  for  United 
States  coal,  Norfolk  may  lose  some  of  this  competitive  advantage.  There  are 
currently  several  vessels  calling  on  Norfolk  that  must  carry  partial  loads  to 
navigate  through  the  existing  channels.  Since  the  majority  of  the  cargo  pass¬ 
ing  through  Norfolk  is  high  in  volume  and  low  in  price,  the  efficient  use  of 
shipping  is  crucial  to  bring  profits.  Unless  the  harbor  is  deepened,  future 
deep-draft  vessels  may  be  forced  to  use  other  ports. 


Proposed  Channel  Improvements 

5-  The  proposed  improvements  to  channels  and  anchorages  approaching 
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Norfolk  Harbor  are  shown  in  Figure  2  and  described  as  follows  (bSAED, 
1980) : 


Increasing  the  depth 
below  mean  low  water 


of  Thimble  Shoal  Channel  from  45  to  55  ft 
over  its  existing  1,000-ft  width. 


b. 


Increasing  the  depth 
below  mean  low  water 
the  coal  terminal  at 


of  Norfolk  Harbor  Channel  from  45  to  55  ft 
over  its  existing  800-  to  1,500-ft  width  to 
Lamberts  Point . 


c.  Increasing  the  depth  of  the  channel  to  Newport  News  from  45  to 
55  ft  below  mean  low  water  over  its  existing  800-ft  width  to  the 
coal  terminal  at  Newport  News. 


d.  Dredging  a  new  channel,  referred  to  as  the  Atlantic  Ocean  Chan¬ 

nel,  off  Virginia  Beach  to  a  depth  of  5?  ft  1  ov  mean  low  water 
and  a  width  of  1,000  ft  over  a  length  of  10.-  lies. 

e.  Constructing  three  fixed-mooring  anchorage  f  'ities,  each  ca¬ 
pable  of  accommodating  two  large  vessels  sim  aneously. 


f.  Increasing  the  depth  of  the  Elizabeth  River  he  Southern 

Branch  of  the  Elizabeth  River  between  Lamber  int  (river 

mile  9)  and  the  Norfolk  and  Western  Railway  bridge  (river 
mile  15)  from  40  to  45  ft  below  mean  low  water  over  its  existing 
375-  to  750-ft  width . 


g.  Increasing  the  depth  of  the  Southern  Branch  of  the  Elizabeth 
River  between  the  Norfolk  and  Western  Railway  Bridge  (river 
mile  15)  anil  the  U.  S.  Routes  460  and  13  highway  crossing  (river 
mile  17.5)  from  35  to  40  ft  below  mean  low  water  over  its  exist¬ 
ing  250-  to  500-ft  width,  and  providing  a  new  800-ft  turning 
basin  at  the  terminus  of  the  channel  improvement. 


6.  The  depths  listed  above  are  project  depths  and  do  not  include  al¬ 
lowable  overdepth  dredging.  The  actual  depths  for  the  proposed  new  channels 
with  the  required  overdepth  dredging  should  be  3  ft  deeper.  Previous  deepen¬ 
ing  projects  in  the  vicinity  also  had  provisions  for  overdepth  dredging.  The 
model  study  of  the  proposed  channel  deepening  used  the  existing  channel  depths 
as  determined  by  the  most  recent  surveys  for  the  base  condition  with  the  plan 
condition  using  the  project  depths  plus  the  3-ft  overdepth  dredging. 


Purpose 

7.  The  purpose  of  this  study  was  to  investigate  the  impact  of  the  pro¬ 
posed  channel  deepening  on  the  hydrodynamic  characteristics  of  Chesapeake  Bay 
The  study  was  designed  specifically  to  determine  what  changes  in  tidal  eleva¬ 
tions,  current  velocities,  and  salinities  could  be  attributed  to  the  proposed 
channel  deepening. 
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Scope  of  Tesli 


8.  The  hydraulic  model  study  consisted  of  two  parts.  The  first  was  a 
series  of  four  steady-state  tests  (constant  discharge  and  cosine  tides)  de¬ 
signed  to  study  base  versus  plan  differences  in  tides  and  current  velocities. 
Both  base  and  plan  geometries  were  tested  under  medium  and  high  tide  range  and 
freshwater  discharge  conditions.  The  boundary  conditions  and  sampling  pro¬ 
cedures  for  the  steady-state  tests  were  dictated  by  the  needs  of  numerical 
models  at  the  U.  S.  Army  Engineer  Waterways  Experiment  Station  (WES)  for  sub¬ 
sequent  studies  of  sediment  transport  and  shoaling  in  the  dredged  channels  and 
neighboring  bottom  areas. 

9.  The  second  part  of  the  study  was  a  dynamic  test  (variable  discharge 
and  variable  tides)  designed  to  predict  base  versus  plan  differences  in  salin¬ 
ity  response.  A  2-1/2-year  variable  hydrograph  was  used  with  a  repetitive 
28-lunar-day  variable  tide  for  both  base  and  plan  geometries.  The  ocean 
source  salinity  was  the  same  for  both  steady-state  and  dynamic  tests. 
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PART  1 1  :  CHESAPEAKE  BAY  MODEL 


Physica  1 


Model  Desci 


10.  The  physical  model  ol  Chesapeake  Bay  is  located  on  Kent  Island  in 
Matapeake,  Maryland.  The  model  is  an  8.6-acre  fixed-bed  model  molded  m  iun- 
crete  to  conform  to  the  most  recent  National  Ocean  Survey  charts.  At  the  time 
of  this  study,  all  major  ship  channels  had  been  molded  with  the  proposed  SO-lt 
channels  leading  into  Baltimore  and  the  existing  channels  elsewhere.  Channels 
in  the  James  and  Elizabeth  Rivers  and  Thimble  Shoal  area  of  the  Lower  Bay  were 
molded  to  correspond  to  prototype  information  collected  as  late  as  1981.  The 
molded  area  of  the  model  extends  from  approximately  JO  miles  offshore  in  the 
Atlantic  Ocean  to  the  heads  of  tide  for  all  tributaries  emptying  into  the 
Chesapeake.  The  entire  length  of  the  Chesapeake  and  Delaware  (C6D)  Canal  and 
a  portion  of  Delaware  Bay  are  also  molded.  Overbank  geometry  is  reproduced  tc 
the  +20  ft  contour.  Model  limits  are  shown  in  Figure  3. 

11.  The  hydraulic  model  was  designed  based  on  the  equality  of  Froude 
numbers,  model  to  prototype,  reflecting  similitude  of  gravitational  effects. 
Geometric  scales  of  the  model  are  1:1,000  horizontally  and  1:100  vertically, 
reflecting  a  distortion  ratio  of  10:1.  These  dimensions  and  Froudian  model 
laws  defined  the  following  mode l -to-prototype  ratios: 


Cha racter i st  i c 
Vertical  length 
Horizontal  length 
Slope 
T  ime 

Velocity 

Volume 

Discharge 


Rat  i  o 
1:100 
1:1,000 
10:  1 
1  :  100 
1  :  10 

1 : 100,000,000 
1 : 1 ,000,000 


The  model-to-prototype  ratio  for  salinity  is  1:1.  This  is  the  general  prac¬ 
tice  for  distorted-scale  models. 

12.  The  model  was  designed  and  equipped  so  that  selected  prototype 
boundary  conditions  could  be  simulated  and  the  model  response  to  these  condi 
tions  recorded.  A  discussion  of  appurtenances  necessary  to  generate  and  re¬ 
cord  test  boundary  conditions  follows. 


gure  3.  Model  limits  and  automatic  water-level  detector  locations  with 
their  computer  number  designations 


Comp  lit  c  r  Fa  r  i  1  i  tics 

13.  The  Chesapeake  Bay  model  is  equipped  with  live  minicomputers  that 
perform  a  variety  of  tasks  ranging  from  model  control  and  data  logging  to  the 
analysis  and  graphical  display  of  model  test  data.  These  include  a  Texas  In¬ 
struments  (Tl)  960,  a  T1  980,  an  International  Business  Machines  (IBM)  9110,  a 
Digital  Equipment  Corporation  (DEC)  PDF  11/23,  and  a  DEC  PDF  11/94. 

14.  The  Tl  960  is  a  64K  minicomputer  used  solely  for  model  control  and 
data  acquisition.  It  is  equipped  with  a  2.5-inegabyte  magnetic  disc  that  con¬ 
tains  all  necessary  system  software  to  compile  and  run  the  model  control  com¬ 
puter  program.  It  also  is  equipped  with  two  250K  flexible  disc  drives  that 
are  used  to  receive  data  from  75  different  flowmeters  and  water-level  de¬ 
tectors  throughout  the  model.  Output  from  these  devices  is  displayed  on  a 
cathode-ray  tube  (CRT)  or  hardcopy  terminal  where  it  can  be  observed  at  the 
same  time  that  it  is  recorded  on  flexible  disc.  Through  the  model  control 
terminal,  an  operator  can  interactively  observe  model  operations  by  displaying 
values  from  any  combination  of  model  control  devices. 

15.  The  Tl  980  is  a  56K  minicomputer  used  primarily  tor  data  analysis. 
It  lus  the  same  access  to  the  magnetic  disc  and  flexible  discs  as  does  the 

TI  960.  In  addition,  it  can  interface  with  a  300  card-per-minute  card  readei , 
a  9-track,  800-bpi  (bytes  per  inch)  magnetic  tape  drive,  and  a  Versatec  elec¬ 
trostatic  printer/plotter.  Graphics  for  this  report  were  partially  supplied 
by  the  Versatec  machine. 

16.  The  IBM  5110  is  a  64K  minicomputer  that  has  access  to  twin  250K 
flexible  disc  drives  and  Tektronix  plotting  peripherals.  It  uses  an  APL  key¬ 
board  and  is  used  in  various  data  editing  and  analysis  tasks. 

17.  Data  logging  at  the  model  can  be  accomplished  through  the  use  of 
the  PDP  11/23  minicomputer.  The  PDP  11/23  has  64K  of  main  memory  and  can 
store  5.8  megabytes  of  data  on  a  Winchester  magnetic  disc.  Although  it  is 
being  developed  as  a  self-contained  data  logging  system,  it  is  also  used  for  a 
variety  of  small  data  analysis  tasks. 

18.  The  largest  of  the  model  minicomputers  is  the  PDP  11/44  which  has  a 
full  megabyte  of  main  memory  and  can  service  eight  users  simultaneously.  Pro¬ 
gram  and  data  storage  in  the  system  is  supplied  by  twin  10.4  megabyte  remov¬ 
able  magnetic  discs,  a  1,600-bpi  tape  drive,  twin  250K  flexible  disc  drives, 
and  a  600  card-per-minute  card  reader.  The  PDF  11/44  is  connected  to  a 
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Tektronix  graphics  system  that  includes  a  4112  video  graphics  screen,  a  4662 
flat  bed  plotter,  and  a  4612  hard  copy  unit.  The  graphics  in  this  report  were 
largely  supplied  by  the  PDP  11/44  and  Tektronix  system. 

Tide  Generation  and  Measurement 

19.  Source  tides  in  the  model  can  be  generated  by  using  the  primary 
tide  generator  in  the  model  ocean  and  by  using  the  secondary  tide  generator  in 
the  Delaware  Bay  at  the  eastern  end  of  the  C&D  Canal.  Both  tide  generators 
can  be  operated  to  generate  a  repetitive  tide,  or  by  using  the  TI  960  model 
control  computer  a  repetitive  or  variable  tide  can  be  generated. 

20.  The  TI  960  computer  controls  the  source  tides  by  providing  a  con¬ 
tinuously  changing  programmed  voltage  to  the  tide-generating  mechanisms. 

These  mechanisms  consist  of  a  feedback  system  that  is  entirely  self-contained 
and  is  not  dependent  on  computer  feedback  for  adjustment.  The  system  consists 
of  a  tide  control  amplifier  that  conditions  the  computer  signal  and  a  bubble 
tube  positioner  that  senses  the  water-level  position  and  positions  the  hydrau¬ 
lically  controlled  inlet  and  outlet  gates  via  a  hydraulic  pilot  amplifier. 
Figure  4  is  a  schematic  of  the  tide-generating  system.  A  more  detailed  de¬ 
scription  of  the  tide  generators  can  be  found  in  Scheffner  et  al.  (1981). 

21.  Water-surface  elevations  throughout  the  model  can  be  measured  both 
manually  by  75  distributed  point  gages,  and  automatically  by  22  water-level 
detectors  (WLD) ,  which  report  their  individual  water  levels  to  the  computer 
where  they  are  stored.  Water-level  elevations  are  monitored  both  manually  at 
the  ocean  and  automatically  throughout  the  model  during  testing  (Figure  3). 
Manual  measurements  at  the  ocean  were  used  to  check  automatic  devices. 

Freshwater  Inflow 

22.  The  Chesapeake  Bay  model  is  capable  of  reproducing  a  variable  hydro¬ 
graph  freshwater  inflow  through  the  use  of  positive  feedback  control  of  river 
discharges.  Fresh  water  normally  enters  the  model  at  21  independent  inflow 
points  representing  the  major  tributaries  of  the  prototype.  During  the  Nor¬ 
folk  test,  however,  an  additional  inflow  was  needed  on  the  Elizabeth  River  to 
more  accurately  simulate  the  inflows  in  the  lower  James  River  watershed.  Or¬ 
dinarily,  Elizabeth  River  flows  enter  the  model  through  the  inflow  on  the 
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THE  WATER  SURFACE  OF  THE  MODEL  (A)  IS  APPROXIMATELY  5  FT  HIGHER  THAN  RETURN  SUMP  (B)  AND 
10  FT  LOWER  THAN  SUPPLY  SUMP  <C)  BECAUSE  OF  THESE  DIFFERENCES  IN  WATER-SURFACE 
ELEVATIONS.  THE  FLOW  OF  WATER  FROM  THE  MODEL  INTO  THE  RETURN  SUMP  AND  OUT  OF  THE  SUPPLY 


Nansemond  River.  Figure  5  is  a  map  of  the  bay  showing  the  positions  of  the 
discharge  points.  The  Susquehanna  River  required  two  inflow  systems  (num¬ 
bers  15  and  22)  due  to  the  range  of  freshwater  inflow.  As  shown  in  Figure  5, 
both  these  systems  lead  to  the  same  discharge  point. 

23.  Although  more  than  100  separate  tributaries  empty  into  Chesapeake 
Bay,  only  22  rivers  were  chosen  to  represent  the  total  combined  bay  discharge. 
Providing  a  separate  discharge  point  at  each  minor  tributary  is  impractical. 
The  sophisticated  plumbing  and  equipment  necessary  for  each  inflow  are  expen¬ 
sive  and  require  specialized  maintenance.  Many  of  the  tributaries  provide  in¬ 
finitesimal  flows,  immeasurable  with  the  present  system  used.  These  flows  are 
summed  with  the  nominal  discharge  of  the  closest  of  the  22  chosen  tributaries 
to  provide  a  representative  and  well-balanced,  as  well  as  a  cost-effective, 
inflow  distribution. 

24.  Flow  is  controlled  at  the  discharge  points  by  an  arrangement  of 
solenoid-controlled  discharge  ports  with  graduated  orifices.  An  electrical 
signal  causes  a  solenoid  to  be  activated,  fully  opening  a  discharge  port.  The 
configuration  of  the  graduated  orifices  is  such  that  4,096  combinations  of 
open  and  closed  ports  provide  a  range  of  flows  which  can  be  stepped  from  the 
smallest  measurable  flow  to  the  individual  tributary's  maximum  flow.  These 
arrangements  of  ports  are  called  digital  valves. 

25.  Flow  from  the  digital  valves  to  the  model  is  monitored  continuously 
by  bearingless  flowmeters  of  varying  ranges  that  can  be  used  alone  or  in  com¬ 
binations  to  provide  accurate  flow  readings  covering  the  full  range  of  a 
tributary's  discharge.  The  flowmeters  use  fiber  optics  to  count  the  revolu¬ 
tions  of  a  water-driven  rotor.  Optical  pulses  are  translated  to  electrical 
pulses  that  are  summed  on  an  arrangement  of  totalizers  and  latches  on  a  coun¬ 
ter  card.  This  counter  is  strobed  at  predetermined  intervals  which  causes  the 
summed  value  of  pulses  to  be  transferred  to  a  transmitter  and  the  latches 
cleared  for  the  next  summation.  During  this  study,  inflows  were  strobed  every 
18  sec,  or  every  half  hour  of  prototype  time;  thus  the  values  available  for 
flow  calculations  are  not  instantaneous  flows  but  are  18-sec  averages. 

26.  Each  flowmeter  in  use  has  its  own  transmitter  where  the  binary  to¬ 
tals  from  the  counter  are  translated  to  ASCII  code  and  transmitted  to  the 

TI  960  model  control  computer  in  serial  form  via  a  hard-wired,  cascading,  mul¬ 
tiplexed  20-ma  current  loop  communication  system. 

27.  In  the  computer,  the  pulse  totals  are  transformed  to  flows  using  a 
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Figure  5 .  Freshwater  inflow  and  sewage  treatment  plant  locations 


linear  regression  derived  f i om  pielest 
calibrations  for  each  f lowmeter/digital 
valve  combination.  The  computer  then 
compares  these  18-sec  averaged  flows 
with  the  desired  flows  for  each  flow¬ 
meter  and  determines  whether  the  digita 
valve  setting  should  be  adjusted  at  any 
of  the  inflows.  The  feedback  system  is 
activated  when  there  is  a  discrepancy 
between  the  desired  flow  and  the  actual 
flow.  Discharge  ports  are  opened  or 
closed  to  adjust  the  actual  flow  toward 
the  desired  flow.  Time-averaged  dis¬ 
charges  controlled  in  this  manner  re¬ 
main  very  close  to  the  desired  hydro- 
graph  step  values.  Some  overshooting 
of  flows  may  occur  at  the  beginning  of 
each  hydrograph  step  as  the  computer 
adjusts,  but  flows  are  generally  stable 
within  just  a  few  update  cycles.  This 


Figure  6.  Typical  inflow  system 


varies,  of  course,  with  the  magnitude  of  the  step  change  as  well  as  with  each 


individual  flowmeter/digital  valve  combination.  Figure  b  shows  a  typical  in¬ 


flow  system  used  in  this  test. 


Sewage  Treatment  Plants 


28.  During  this  study,  changes  in  discharges  ol  the  rivers  were  sup¬ 
plemented  by  modeling  major  sewage  treatment  plants  (STP's)  in  the  Noitolk 
area.  Figure  5  shows  the  locations  of  the  live  STP's  modeled  tor  the  study 
Table  1  is  a  list  of  the  stations  with  their  geographical  locations  Dis¬ 
charges  A  and  B  are  located  on  the  James  River,  discharges  C  and  I)  are  locator 
on  the  Elizabeth  River,  and  discharge  F  is  located  near  Little  Creek  in  the 
Lower  Chesapeake  Bay.  Engineering  drawings  of  the  STP's  were  consulted  and 
the  discharge  points  were  located  as  near  as  possible  to  those  in  the  proto¬ 
type.  Fresh  water  was  used  as  the  discharge  medium  in  all  modeled  STP's. 

Since  ncar-ficld  flow  patterns  are  not  easily  modeled  in  a  distorted-scale 
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model,  and  Mince  a  near-field  flow/ 
density  structure  was  not  an  object  of 
concern,  no  attempt  was  made  to  repro¬ 
duce  tlu*  injection  methods  oi  prototype 
Brass  ditlusers  were  fitted  to  nylon 
tubing  lot  out  tails.  Figure  7  shows  a 
typical  portable  inflow  as  used  m  the 
Sot  I  <>lk  Harbor  study.  A  const  ant-head 
tank,  approx (male  1 y  10  II  above  the 
node!  ,  tec. Is  an  array  of  adjustable  ro¬ 
tameters  whi.ti  in  tuin  teed  t  tie  nylon 
d  i  si  lia  t  ge  1  i  ties  1  1  ow  rates  through 
the  Sir's  were  adjusted  manually,  as  re 
•  pitied,  t  h  r " u  ghe  u  t  each  hydrograph  test 


Sinr\  Nwi  leai  Power  Plant 


.■>».  The 


> 1  l ng  wale r  divers i on 


Figure  7.  Tvpi.  it  SIT 


tr.im  t  tie  Surry  Nucieai  Power  Plant  was 
simulated  in  model  testing  tor  the 


tiisl  tune  in  t  tie  N.ut.lk  H.nhor  study.  1'hc  plant  is  located  at  Hog  Island  on 
t  he  southern  stioie  ot  i  lie  lanus  Hivei  (figuies  t  and  8).  Virginia  Electric 
I’owei  t.  imp.  inv  iVF.Pi't')  .'ttiiials  have  estimated  l  tie  diversion  at  1,740  cts  con¬ 
tinuously  even  w  tie  ii  t  lie  plant  is  not  on  1  i  ne  A  diversion  o  t  this  ma  gn  1 1  ude 
is  several  times  gieatei  than  l  tie  total  James  Kiver  inflow  during  drought  Con¬ 
di  (  ions  so  its  i  in  I  us  i  mi  was  deemed  net  c.ss.ny  for  t  h  i  s  study. 

to  Engineering  dtawings  and  hvd rograph i i  surveys  were  consulted  so 
that  the  intake  and  out  tall  geometries  as  well  as  their  nearshore  approach 
i  flannels  were  modeled  to  closely  approximate  prototype  conditions.  figure  8 
shows  the  intake  and  out  tall  lot  at  ions  near  Hog  Island  and  the  relative  close¬ 
ness  to  the  nearest  salinity  sampling  stations.  During  the  model  test,  the 
1,740-t  Is  (prototype)  lonstant  bypass  was  maintained  by  reading  a  discharge 
rotameter  ami  adjusting  a  throttling  valve  which  was  attached  to  a  constant 
rpm  centrifugal  pump  (Figure  4).  An  identical  calibrated  backup  pump  was 
available  in  case  the  primary  pump  failed  during  testing.  No  significant  dif¬ 
ferences  in  pump  performance  were  not  i  <  e<1  between  the  base  and  plan  tests. 
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.<1.  A  constant  source  salinity  is  provided  to  the  model  <<<  c,in  F,  in.iiii- 
tainiug  the  supply  sump  (Figure  A)  .it  the  drsiied  s.ilinitv.  I'ti  i  :  is  .11- 
comp  l  i  sited  by  adding  suft  n  lent  <|itant  it  irs  o  l  sat  in  at  H  In  i  nr  I  app  i  .  i  m.i  t  •  ,  \ 
280  ppt  )  to  the  return  sump  where  it  is  mixed  t--  tin  desiied  s  1 1 1  .  •  -  ilii.it'. 

prior  to  being  circulated  to  the  supply  sump  and  then.,  t..  the  im  d*  I  .■>  eat.  . 

The  saturated  brine  solution  is  obtained  by  injecting  fresh  watei  into  a  bed 
of  granular  salt  ( N a C 1  I  from  which  it  is  latei  lelea-ed  to  the  let'irn  •  ump  .u 
measured  quantities.  The  saltwater  supply  system  is  i.ip.ihlr  maintaining 
sourer'  salinity  to  within  0.2  ppt  of  that  desired  m  s t ea.h  -  a ! a t e  londiti  us 
and  to  within  0.8  ppt  of  that  desired  during  variable  hvdrographir  r oud  i  i  ions. 


Bubbler  System 


12.  The  bubbler  system  in  the  model  is  designed  to  create  a  more  real¬ 
istic  vertical  salinity  distribution.  Since  nonast ronomi ca 1  mixing  energy 
(primarily  wind)  is  not  easily  modeled  by  tides  and  supplemental  roughness 
alone,  it  was  necessary  to  add  the  bubbler  system  m  order  to  maintain  closer 
agreement  with  the  vertical  salinity  distribution  in  the  prototype. 

33.  The  model  bubbler  system  consists  of  a  network  of  copper  tubing 
placed  along  the  axis  of  the  bay  and  its  major  tributaries  (Figure  101.  The 
tubing  is  charged  with  a  constant  air  pressure  and  releases  bubbles  into  the 
water  column  at  a  constant  prescribed  flow  rate.  Throughout  any  given  test 
the  bubble  flow  rate,  air  pressure,  and  bubbler  depth  in  the  water  are  moni¬ 
tored  for  consistency. 


Current  Meters 

34.  In  this  study,  both  magnitude  and  direction  of  currents  were  sam¬ 
pled.  Current  magnitude  measurements  were  made  with  miniature  Price-type 
meters  (Figure  11).  The  center  line  of  the  model  cups  on  the  meter  was  about 
0.04  ft  above  the  bottom  of  the  meter  frame.  The  overall  width  of  the  meter 
was  about  0.1  ft  in  the  model,  representing  a  horizontal  width  of  about  100  ft 
in  the  prototype.  Therefore  distortion  of  the  horizontal  area  (model  to  proto 
type)  resulted  in  model  velocities  averaged  over  a  much  larger  area  than  those 
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Figure  11.  Miniature  Price-type 
velocity  current  meter 


of  the  prototype  point  observations.  The  same  was  true  for  the  vertical  area 
since  the  height  of  the  cups  on  the  meter  was  equivalent  to  about  4.0  ft  pro¬ 
totype.  Current  speeds  were  obtained  by  counting  the  number  of  revolutions 
the  meter  made  in  a  10-sec  interval  which  was  equivalent  to  about  17  min  in 
the  prototype.  The  meters  were  calibrated  frequently  to  ensure  the  ai.iur.uv 
of  measurements  and  were  capable  of  measuring  actual  speeds  as  low  as  about 
0.0:1  lps  (0.  1  fps  prototype).  Accuracy  of  these  meters  was  about  ±0 .  1  5  i  ps 
( p  rot o t  ype  )  . 


15.  Current  illicit  unis  were  recorded  hourly  t  o  the  nearest  10  deg  by 
using  a  direction  vane  and  direction  rose  as  shown  in  Figure  1C.  As  was  true 
with  the  Price-type  meters,  the  size  of  the  vane  results  in  direction  readings 
that  are  an  average  value  over  its  2-ft  (prototype)  height.  In  most  cases, 
current  direction  changes  are  tairly  slow  so  technicians  could  maintain  read¬ 
ing  precision  to  within  10  deg. 
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Figure  12.  Current  direction  vane  and  compass  rose 
Vacuum  S  amp  Id  n  g  Sy stem 


36.  Salinity  samples  for  this  test  were  taken  using  a  vacuum  aspiration 
system.  A  series  of  vacuum  pumps  provide  a  continuous  vacuum  to  three  valve 
manifolds.  A  total  of  16  valves  control  a  varying  number  of  cotidal  stations 
arranged  such  that  sampling  times  are  never  more  than  1/2  hr  (prototype)  from 
the  desired  slack-water  sampling  time.  In  each  case,  the  sample  is  pulled  at 
the  same  time  relative  to  slack  water.  The  sampling  station  consists  of  a 
sampling  probe  built  from  a  number  of  separate  copper  tubes  soldered  together 
to  form  a  multidepth  probe.  These  copper  tubes  are  attached  to  short  lengths 
of  plastic  tubing  that  lead  to  individual  10-ml  test  tubes.  Vacuum  is  pro¬ 
vided  to  the  tubes  by  a  vacuum/overflow  jar  that  provides  an  even  vacuum  to 
all  lines.  This  system  has  proven  effective  in  taking  a  very  large  number  of 
samples.  The  filled  test  tubes  are  removed  manually  from  the  sampling  device 
between  sampling  times  and  placed  in  special  racks  for  later  analysis. 


Salinity  Testing  System 

37.  Beckman  RA-5  solumeters  are  the  primary  instruments  for  making  con¬ 
ductivity  measurements  at  the  model  (Figure  13).  These  meters  use  a  salinity 
probe  shaped  like  an  eyedropper  into  which  a  sample  is  drawn  and  the  salinity 
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Figure  13.  Beckman  solumeter 


Figure  14.  Salinity  data  logger 
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is  read  from  an  analog  meter.  Analog  meters  provide  much  opport un i tv  lor 
error  in  reading  and  transcription,  and  ev<  tlie  most  experienced  operators 
are  unable  to  be  more  accurate  than  2  percent  ol  full  scale,  which  is  0.8  ppt 
The  solumeters  provide  a  voltage'  output  of  0  to  100  mv  which  is  proportional 
to  0  to  100  percent  of  full  scale.  Tins  voltage  is  used  t<>  drive  a  digital 
voltmeter  that  measures  to  within  ±0.1  mv  giving  a  resolution  of  0.1  percent 
of  full  scale  or  0.04  opt.  In  calibrating  the  meters,  it  was  found  that  a 
given  salinity  standard  could  be  repeated  to  within  approximately  0 .  .3  percent 
or  0.12  ppt.  This  is  in  an  ideal  situation  where  temperature's,  conductivity, 
and  probe  residence  time  are  carefully  control  Ic'd.  Other  tests  performed  on 
the  meters  indicate  that  most  samples  can  be  relied  upon  to  within  approxi¬ 
mately  0 .  S  ppt  .  There  are  some  cases  where  operator  emit  or  electronic  proh 
1  ems  can  cause'  larger  deviations,  but  these  are  detectable  and  can  be'  iso¬ 
lated  from  the-  data  set  or  corrected. 

.38.  The'  data  logging  system  (Figure*  141  used  tor  tins  study  outers  the 
values,  converted  hv  the'  digital  voltmeters,  on  cassette  tape  in  ASCII  code 
where  it  can  he'  processed  by  the  11  ‘180  minicomputer  tor  storage.  With  each 
sample  value,  pertinent  information  such  is  depth,  time,  tide,  and  station 
name  is  added  to  the*  record.  This  can  all  no  accomplished  with  a  minimum  ef¬ 
fort  on  the  part  of  the  meter  operator.  Direct  entry  of  values  on  cassette 
tape  precludes  the-  need  for  keypunch  i  tig  and  the  possibility  of  m  l  s  i  nt  e  rp  re  t  u- 
liott  with  each  transfer  to  a  different  medium. 


PART  III.  STEADY- ST ATL  TIDE  AND  VELOCITY  TESTS 


Test  Conditions 


39.  The  purpose  of  the  steady-state  tests  in  the  Norfolk  Harbor  study 
was  to  examine  the  impact  of  channel  deepening  on  tidal  elevations  and  current 
velocities  in  the  lower  Chesapeake  Bay  and  James  River  regions.  The  data  col¬ 
lected  from  these  hydraulic  model  tests  were  to  he  used  in  numerical  hydro- 
dynamic  models  at  WES  under  a  separate  agreement  with  the  l).  S.  Army  Engineer 
District,  Norfolk,  to  study  sediment  transport  and  shoaling  characteristics  of 
the  dredged  Norfolk  shipping  channels  aim  adjacent  estuary  bed  areas.  The 
proposed  improvements  to  Norfolk  Harbor  and  its  approach  channels  are  de¬ 
scribed  in  PART  I  and  are  shown  in  Figure  2. 

40.  The  term  "steady  state"  indicates  that  the  hydraulic  model  was  op¬ 
erated  under  conditions  of  constant  river  inflow  and  a  repetitive  cosine  tide 
at  the  model  ocean.  The  testing  scheme  consisted  of  a  base  condition  in  which 
the  model  bathymetry  simulated  present  prototype  conditions  with  the  Norfolk 
channels  molded  to  depths  reported  from  a  l' .  S.  Army  Corps  of  Engineers  hydro- 
graphic  survey  in  1980,  and  a  plan  condition  in  which  the  channels  were 
deepened  to  their  proposed  project  depths  plus  3  ft  for  overdepth  dredging. 

The  tide,  velocity,''"  and  direction  testing  for  both  base  and  plan  conditions 
was  performed  under  the  following  four  test  scenarios: 

Test  I:  The  model  was  operated  under  a  repetitive  cosine  tide  vary¬ 
ing  from  +2.40  to  -2.40  ft  N0YP'"v  at  the  Atlantic  Ocean 
control  station.  This  approximates  a  spring  tide  condition 
in  Chesapeake  Bay.  The  total  bay  freshwater  river  inflow 
was  a  constant  200,000  cfs  which  represents  a  relatively 
high  flooding  condition. 

Test  2:  The  total  bay  discharge  remained  at  200,000  els;  however, 
the  repetitive  cosine  tide  varied  from  +1.50  to  -1.50  ft 
NGVD  which  approximates  a  neap  tide  condition. 

Test  3:  The  total  hay  discharge  was  70,000  cfs  which  simulates  the 
long-term  average  flow  into  the  Chesapeake  Bay  from  all  of 

Velocity  can  be  rigorously  defined  as  a  vector  quantity  containing  speed 
and  direction.  It  is,  however,  commonly  used  as  a  scalar  quantity  describ¬ 
ing  current  speed  in  either  a  flood  or  ebb  direction.  Hereafter,  velocity 
will  be  used  as  a  description  of  current  speed  unless  otherwise  noted. 

**  All  elevations  cited  herein  are  in  feet  referred  to  the  National  Geodetic 
Vertical  Datum  (NGVD),  though  on  some  figures  Mid  plates  "ft  ms  I "  is  used. 
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+  2.40  to  -2.40  ft  was  generated  at  the  ocean. 


Test  4:  The  total  hay  discharge  remained  at  70,000  its  and  the  ap¬ 
proximate  neap  tide  varying  from  ♦ 1 . SO  to  -1.S0  (t  was 
generated  at  the  ocean. 

41.  Steady-state  discharges  tor  each  ot  the  21  model  river  mllows  are 
summarized  in  Table  2.  In  addition  to  the  river  inflows,  live  sewage  treat¬ 
ment  plants  and  the  Surry  Nuclear  Power  Plant  cooling  water  diversion  were 
operated  throughout  the  four  steady-state  tests.  Locations  of  the  sewage 
treatment  plants  and  the  Surry  plant  are  shown  in  figure  1  with  their  loordi- 
lutes  and  si eady- si  at e  disiharge  rates  given  in  Table  1. 

42.  During  both  the  s t eady-s t a t e  and  dynamic  tests,  tides  were  gener¬ 
ated  only  at  the  Atlantic  Ocean  boundary.  Tile  model  is  capable  of  generating 
tides  at  the  Delaware  Ray  end  of  the  Oil)  Lana  I  hut  this  was  not  done  dining 
the  Norfolk  study.  1  he  Delaware  Ray  tide  source  is  sufficiently'  distant  from 
the  study  area  to  make  its  inclusion  insignificant  and  certainly  not  cusl- 

e I  t  ect  i ve . 

41.  Ihe  source  salinity  lo  the  entire  model  study  was  12 . h  ppt  .  Con¬ 
trol  o t  the  source  salinity  -tilling  the  steady-slate  tests  was  considered  good 
with  any  minor  variations  being  incapable  of  causing  any  observable  base 
versus  plan  differences  in  tides  or  velocities.  A  complete  document  at i on  ot 
all  model  boundary  conditions  during  these  tests  is  available  (Waltz  and 
Morton  1482). 


fide  Test  Data  and  Results 

44.  During  steady-state  testing,  model  tide  observations  were  made  at 
the  lb  stations  shown  in  figure  IS  and  at  one  in  the  ocean.  The  17  stations 
are  listed  in  Table  I.  fide  elevations  were  recorded  hourly  (every  Tb  sec 
real-time)  over  two  tide  cycles  or  24.84  prototype  equivalent  hours.  These 
tide  measurements  were  made  both  manually  using  graduated  point  gages  attached 
to  portable  tripod  supports  and  automatically  using  electronic  water-level 
detectors.  Data  from  the  water-level  detectors  were  stored  on  flexible 
diskettes  via  the  model  control  computer. 

45 .  Tide  station  locations  were  determined  by  the  requirements  of  the 
numerical  hydrodynamic  models  to  be  used  in  the  sediment  transport  studies  at 
WES.  Major  boundaries  of  the  two  numerical  hydrodynamic  models  were  at  range 
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CBOO  on  the  cdbt,  CB01  on  the  north,  range  JG01  at  the  mouth  of  the  James 
River,  and  range  JG03  at  the  upstream  boundary  on  the  James.  To  adequately 
describe  cross-sectional  water-surface  slopes,  three  tide  stations  were  lo¬ 
cated  across  range  CBOO  at  the  mouth  of  Chesapeake  Bay,  and  three  stations 
were  positioned  across  range  CB01.  Since  the  James  River  boundary  ranges  had 
relatively  short  widths,  their  cross-sectional  water-surface  slopes  were  es¬ 
sentially  horizontal  in  the  model  and  could  be  described  by  a  single  tide 
station. 

46.  During  the  base  and  plan  steady-state  tests,  over  3,400  discrete 
hourly  tide  elevation  data  points  were  collected  from  the  17  manual  sampling 
stations.  The  data  set  is  complete  with  the  exception  of  sta  TS0005  which  is 
missing  from  base  tests  2  and  3.  From  each  station  for  each  test,  23  con¬ 
secutive  hourly  tide  elevations  constituting  two  complete  tidal  cycles  were 
fitted  to  the  following  expression: 


where 


hit  1 


A  +  a 
o 


cos 


2mat 

360 


2n0\ 

360  / 


(1) 


hft)  =  tide  height,  ft 
t  =  t ime ,  hr 

=  mean  water  level,  ft,  of  hourly  tides 
a  =  amplitude,  ft 

uj  =  M(  constituent  angular  velocity,  23.984104  deg/hr 
0  =  phase  angle  or  epoch,  deg,  from  zero  lunar  hour 

47.  Resulting  cos ine- f i t ted  curves  for  the  observed  tides  are  shown  in 
Plates  1-17,  and  the  computed  mean  water  level  (A  ),  amplitude  (a),  and  phase 
angle  (0)  are  listed  in  Tables  4A-4D  along  with  base  to  plan  comparisons.  The 
value,  R- v'“ 2  ,  shown  in  Plates  1-112  is  a  statistic  that  measures  the  goodness 
of  fit  of  the  cosine  curve  through  the  individual  data  points  with  a  value  of 
1.0  being  an  exact  fit. 

48.  Manual  point  gages  are  graduated  to  0.001  ft  (0.1  ft  prototype 
equivalent)  and  can  be  estimated  to  the  nearest  0.000.3  ft  (0.03  ft  prototype). 
Thus  it  is  reasonable  to  assume  that  the  accuracy  of  the  manual  point  gages 
for  a  single  tide  measurement  is  ±0.10  ft  (prototype)  depending  upon  a  number 
of  variables  including  the  keenness  of  the  technician's  eye,  the  delicacy  of 
his  or  her  touch  in  lowering  the  point  gage  to  the  water  surface,  and  the  pres 
cncc  or  absence  of  small  water-surface  ripples  at  the  instant  of  measurement. 
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A.  v  ui.u  )  vj  I  tin-  i  wmput  ed  mean  water  !<vrl  and  tidal  amplitude  1 there!  art* 
also  assumed  to  be  10. 10  tt  (prototype). 

44.  "Zero- luiur-hour"  is  the  time  reference  used  to  ensure  that  mea¬ 
surements  made  in  the  model  can  be  compared  with  those  taken  in  the  prototype. 
Zero- l unat -hour  occurs  when  the  moon  is  directly  over  Chesapeake  Bay.  For  the 
model  reveri  I  ic.it  ion  tests  conducted  in  1481,  the  model  clock  in  the  ocean 
tide  control  hut  was  set  so  that  zero- lunar-hour  occurred  at  hour  6.6  (141  deg 
epoch).  This  same  clock  setting  was  also  used  tor  both  base  and  plan  tests  of 
the  Norfolk  Harbor  study.  The  procedure  tor  setting  the  model  clock  was  some¬ 
what  subjective.  A  technician  would  take  a  succession  of  readings  at  the 
ocean  point  gage,  the  time  between  readings  normally  being  between  6  and 
10  sec  l real-time).  Upon  noticing  that  the  peak  of  the  tidal  wave  had  oc¬ 
curred,  the  technician  would  then  note  the  time  on  the  model  clock  and  adjust 
it  to  houi  6 . 6  if  necessary.  The  time  lag  from  the  instant  of  actual  tidal 
peak  occurrence  to  the  instant  that  the  technician  perceives  this  peak  could 
be  as  great  as  17  sec  (14  deg  epoch).  In  general,  no  adjustments  were  made  to 
the  clock  unless  high-water  occurrence  varied  by  0.5  hr  or  more  from  the  de- 
s i red  hou  r  6.0. 

ri0.  Technicians  are  cued  by  a  bell  timing  system  which  rings  for  a  dura 
lion  of  10  sec  at  the  start  of  each  prototype  hour,  and  they  are  instructed  to 
make  their  tide  measurements  as  near  as  possible  to  the  beginning  of  the  bell. 
It  is  believed  that  discrepancies  in  time  of  measurement  could  cause  as  much 
as  an  8-deg  error  in  epoch  when  considering  a  discrete  tide  observation. 

Hence  it  is  reasonable  to  assume  that  the  raw  model  phase  angles  are  accurate 
to  within  i22  deg. 

51.  P 1 a»-mi mis -base  amplitude  differences  for  the  four  tests  are  summa¬ 
rized  in  Table  5.  Very  little  change  in  amplitude  was  noticed  between  base 
and  plan  tests.  From  the  differences  listed  in  Table  5,  one  may  be  able  to 
deduce  trends  indicating  a  slight  increase  in  amplitude  for  the  plan  condition 
during  Test  1,  slight  decreases  in  plan  amplitudes  during  Tests  2  and  3,  and 
essentially  no  change  during  Test  4.  However,  the  largest  change  in  amplitude 
measured  in  the  model  was  -0.08  ft  at  sta  CB0008  during  Test  ).  This  is 
within  the  accuracy  limitations  of  the  model;  thus  any  differences  listed  in 
Table  5  are  considered  due  to  random  measurement  error  inherent  in  the  model 
and  not  a  result  of  channel  deepening. 

52.  Plan-minus-base  mean  water-level  differences  for  all  four  tests  are 
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summarized  in  Table  6.  Computed  mean  watei-level  diffeieiu.es  f  i  uili  the  manual 
point  gages  were  nearly  un i nte rpre tab le  due  to  two  problems  that  existed  in 
the  model  during  testing: 

a.  The  portable  tripods  on  which  the  manual  point  gages  were 
mounted  were  not  stable  enough  for  the  degree  of  accuracy  re¬ 
quired  to  interpret  base-to-plan  mean  water-level  changes.  The 
procedure  for  referencing  these  point  gages  to  a  datum  involved 
lowering  the  point  of  the  gage  to  a  submerged  fixed  pin  whose 
elevation  was  known.  However,  the  pressure  exerted  in  touching 
the  point  to  the  pin  was  apparently  great  enough  in  many  cases 
to  cause  the  entire  tripod  to  tilt  by  as  much  as  0.30  ft  (proto¬ 
type).  This  may  account  for  the  anomalous  values  in  Table  6. 

b.  In  the  ocean  tide  hut,  the  model  is  controlled  by  a  bubble  tube 
positioner  (BTP)  and  a  fixed-point  gage  housed  inside  a  baffle 
pit.  The  main  purpose  of  the  baffle  pit  was  to  prevent  surface 
waves  from  interfering  with  the  tide  sensing  equipment.  Prior 
to  the  start  of  the  Norfolk  Harbor  study,  this  baffle  pit  was 
reconstructed  using  solid  barriers  to  prevent  the  thin  fresh¬ 
water  lens,  which  accumulates  in  the  ocean  as  a  result  of  river 
inflows,  from  being  sensed  by  the  BTP.  This  baffle  pit,  there¬ 
fore,  prevented  the  horizontal  exchange  of  ocean  water  on  the 
outside  with  essentially  stagnant  water  on  the  inside.  Near 
the  completion  of  the  plan  test,  an  experiment  was  performed  in 
which  it  was  discovered  that  the  water  inside  the  baflle  pit 
was  A. 5  ppt  less  saline  than  an  equivalent  column  of  ocean 
water  on  the  outside. 

53.  The  model  computer  controlled  the  level  in  the  model  so  that  the 
mean  water  level  at  the  ocean  station  within  the  baffle  pit  was  0.00  ft.  How¬ 
ever,  the  laws  of  fluid  statics  state  that  existence  of  a  4.5-ppt  density  dif¬ 
ference  will  cause  the  plane  of  the  model  outside  the  baffle  pit  to  be  0.14  ft 
less  than  the  plane  within  the  baffle  pit.  The  only  fixed-point  ,age  other 
than  the  ocean  station  that  was  monitored  during  the  Norfolk  Harbor  tests  was 
sta  3  at  Old  Point  Comfort  on  the  James  River.  It  is  evident  in  Table  6  that 
the  plan-minus-base  mean  water-level  difference  is  -0.14  ft  at  sta  3  just  as 
the  aforementioned  experiment  concluded.  No  data  are  available  from  the  base 
tests  to  confirm  whether  a  density  difference  was  present  in  the  ocean  baffle 
pit;  however,  based  on  the  data  from  sta  3,  it  is  assumed  that  no  such  differ¬ 
ence  existed  at  that  time.  Therefore  it  is  reasonable  to  treat  this  as  a  sys¬ 
tematic  error  and  apply  an  approximate  additive  correction  of  +0.14  ft  to  all 
base-to-plan  mean  water-level  differences  (see  adjusted  mean  water-level  dif¬ 
ferences  in  Table  6). 

54.  Although  the  unstable  tripods  have  created  an  error  which  cannot  be 
filtered  out  of  the  mean  water-level  data,  it  is  interesting  to  note  that  the 
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oveiaii  average  hi  ad  j  us  L  ed  mean  walei-leVel  uilleieiu.es  in  Table  6  was  0.00  It. 
tor  Test  1,  +0.06  tt  tor  Test  2,  +0.01  ft  for  Test  3,  and  -0.05  tt  for  Test  9. 
The  overall  mean  of  all  tour  tests  combined  showed  a  difference  of  0.00  ft. 
Based  on  these  averages,  a  general  inference  tan  he  made  that  on  t  lie  whole, 
channel  deepening  caused  no  changes  in  mean  water  levels  in  this  model  study. 

55.  All  phase  angles  listed  in  Tables  4A-4D  have  been  normalized  to  a 
standard  ocean  station  phase  of  191  deg  and  are  summarized  in  Table  7  along 
with  plan-minus-base  differences.  By  normalizing  all  phase  angles  to  the 
ocean  station,  the  error  in  model  control  of  the  occurrence  of  high  water  at 
hour  6.6  has  been  filtered  out  of  all  measurements.  Any  remaining  basc-to- 
plan  discrepancies  in  phase  can  now  be  attributed  either  to  the  time-of- 
measurement  random  sampling  error,  or  to  the  effects  of  channel  deepening,  or 
a  combination  of  both.  Plan-minus-base  differences  in  Table  7  are  all  within 
the  ±18  deg  t ime-of-measurement  error  hand.  Thus  the  model  evidence  indicates 
that  all  phase  differences  are  due  to  random  sampling  errors  rather  than  the 
result  of  deepening  Norfolk  Harbor  and  its  approach  channels. 

Velocity  Tests 

56.  Current  velocities  in  the  Chesapeake  Bay  model  are  largely  the  re¬ 
sult  of  three  factors:  tides,  freshwater  discharges,  and  salinity  density  gra¬ 
dients.  Of  these,  the  most  important  factor  affecting  current  velocities  is 
the  amplitude  of  the  tide  which  varies  throughout  the  model  depending  on  the 
local  bathymetry.  Freshwater  discharge  plays  a  variable  role  depending  upon 
its  magnitude  and  the  estuary  cross-sectional  geometry,  the  most  significant 
effect  being  in  the  upper  reaches  of  tributaries.  Existence  of  a  salinity  den¬ 
sity  gradient  causes  a  multilayer  flow  system  to  be  formed  in  which  the  denser 
bottom  water  tends  to  move  predominantly  in  the  flood  direction  (upstream! 
while  the  fresher  surface  layer  flows  mainly  in  the  ebb  direction  (downstream). 

57.  Model  velocities  were  measured  at  32  sampling  stations  during  the 
four  steady-state  tests.  Velocity  observations  were  collected  concurrently 
with  the  tide  elevations  discussed  in  the  previous  section.  In  addition,  cur¬ 
rent  direction  data  were  collected  at  26  of  the  velocity  sampling  stations. 

Six  stations  (EH  0203,  0501,  0701,  0901,  EE  0301,  and  WB  0201)  were  located  in 
constricted  channel  areas  where  direction  measurements  were  unnecessary. 
Locations  of  the  velocity  stations  are  shown  in  Figure  16  and  are  listed  in 
Table  8. 
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58.  Station  locations  were  determined  primarily  by  boundaiy  require¬ 
ments  of  the  numerical  hydrodynamic  and  sediment  transport  models  at  WES. 
Velocities  were  sampled  across  ranges  CBOO,  CB01  ,  JG01  ,  and  JG03  which  defined 
transverse  boundaries  in  the  mathematical  models.  Additional  samples  were 
taken  in  the  Norfolk  approach  channels  and  in  the  Elizabeth  River  and  its 
tributaries.  Depending  upon  total  station  depth,  velocities  and  directions 
were  generally  sampled  at  the  surface,  middle,  and  bottom.  At  stations  lo¬ 
cated  within  the  Norfolk  approach  channels,  velocities  were  collected  at  coin¬ 
cident  depths  for  both  base  and  lan  channel  configurations  with  an  additional 
deeper  bottom  depth  sampled  during  the  plan  test  (Figure  17). 


Figure  17.  Typical  velocity  sampling  scheme 


59.  The  testing  equipment  included  miniature  Price-type  velocity  meters 
and  a  small  direction  vane,  both  of  which  are  discussed  in  PART  11.  Velocity 
magnitude  and  direction  measurements  were  recorded  every  hour  (36 
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li.it)  .i  unique  calibration  equation  which  related  the  current  velocity  (proto¬ 
type  fps)  to  the  number  ot  revolutions  traveled  by  the  meter  cups  in  a  10-set 
interval  (real-time).  Each  meter  also  had  a  so-called  " t hresho I d-o t -mot  1  on" 
defined  as  a  critical  velocity  less  than  which  the  meter  would  tease  to  turn. 

At  velocities  less  than  critical,  the  momentum  of  the  moving  water  was  not 
great  enough  to  overcome  the  inertial  and  frictional  forces  ot  the  meter.  The 
t hresho 1 d-of -mot i on  varied  from  meter  to  meter  and  was  generally  0.  10  fps 
(prototype  equivalent)  or  loss. 

t)0 .  Invert  a  inty  in  divide  1  velocity  measurements  results  from  errors 
falling  into  the  following  three  classes: 

a.  Hotmdary  condition  errors--These  are  caused  by  variations  in 
source  salinity,  tide  control,  river  discharge  control,  and 
local  perturbations  in  the  flow  field. 

1).  Instrument  e r rors - -These  result  from  imperfections  in  the  manu¬ 
facture  and  calibration  of  the  velocity  meters.  As  an  example, 
the  mechanical  cups  of  a  meter  may  not  be  centered  symmetrically 
about  its  axis  causing  an  irregular  rotation  even  though  the 
actual  current  may  he  steady  and  uniform. 

c.  Personal  errors- -These  errors  stem  t rom  limitations  ot  the 
human  senses  of  sigiii,  lotnh,  and  hearing,  all  ot  which  are 
used  during  velocity  testing. 

Three  of  the  boundary  condition  errors  (source  salmitv,  tide  control,  and 
river  discharge)  were  maintained  within  acceptable  limits  throughout  the 
steady-state  testing;  and  the  recorded  variations  will  have  a  negligible  ef¬ 
fect  on  ■  reliability  ot  current  velocities.  Local  perturbations  in  the 
I  low  field  are  random  natural  errors  and  tend  to  cancel  out  with  a  large  'lum¬ 
ber  ot  observations.  instrument  errors  were  minimized  hv  screening  all  avail¬ 
able  velocity  meters  ill  the  initial  calibration  process  and  choosing  tor  test¬ 
ing  purposes  only  those  exhibiting  the  highest  quality.  The  calibration  of 
each  meter  was  checked  at  least  once  ,.  '.a;  in  the  laboratory  flume.  The  cali¬ 

bration  equations  are  accurate  to  withi'.  approximately  ±0.15  fps  based  on 
95  percent  confidence  interval  with  47  deorp.  s  of  freedom.  The  most  signifi¬ 
cant  personal  error  introduced  into  the  velocity  measurements  was  in  estimat¬ 
ing  the  tractions  of  a  revolution  that  a  meter  turns.  This  can  fie  reasonably 
estimated  to  the  nearest  0.2  revolutions  or  about  ±0.15  fps.  Thus,  consider¬ 
ing  all  the  sources  of  possible  error,  the  accuracy  of  a  discrete  model  veloc¬ 
ity  measurement  is  approximately  ±0.50  fps  (prototype  equivalent). 
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til.  Ilit*  t  rein  a  i  dimension  oi  nm-  velocity  meter  up  it.  0  .  64  fi  i.u  nidi  i 
in  diameter  and  lias  an  area  of  0.001265  stj  ft.  Due  to  the  distorted  stale  of 
t  lit'  inode  J,  t  h  i  s  corresponds  to  an  otpiivalerit  prototype  area  of  126.5  sq  ft  in 
the  shape  of  an  ellipse  having  a  semimajor  horizontal  axis  of  200  ft  and  a 
seininu  nor  vertical  axis  of  2  ft.  Because  of  l  fi  i  s  scale  effect,  velocities 
measured  in  the  model  cannot  he  interpreted  as  equivalent  prototype  point 
velocities;  instead,  t  iiev  correspond  to  prototype  velocities  averaged  over  the 
area  of  the  equivalent  ellipse. 


Velocity  Data  Analysis 


62.  Over  12,000  discrete  velocity  and  direction  measurements  were  re¬ 
corded  during  the  base  and  plan  tests.  At  each  sampling  depth,  25  consecutive 
hourly  velocity  observations  constituting  two  tidal  cycles  were  fitted  to  the 
following  equation  using  a  harmonic  least  squares  curve  fitting  technique. 


v(t ) 


U  +  u  cos 
o 


/  27twt 
\  360 


2n6\ 

160/ 
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where 

v(t)  =  current  velocity,  fps 
t  =  time,  hr 

=  mean  velocity,  fps  of  the  hourly  observations 
u  =  amplitude,  fps 

w  =  M,(  constituent  angular  velocity,  28.484104  deg/hr 
0  =  phase  angle,  deg 

The  confidence  of  the  computed  velocity  constants  (phase  (0),  amplitude  lu), 
and  mean  (U  ))  is  generally  greater  than  that  of  a  discrete  measurement  since 
all  25  hourly  observations  are  used  in  the  analysis.  The  velocity  constants 
are  relatively  insensitive  to  a  few  outliers. 

65.  Individual  velocity  observations  and  cos ine- f i t t ed  curves  for  both 
the  base  and  plan  conditions  are  shown  in  Plates  18-112.  Velocity  constants 
are  listed  in  Tables  9A-9D.  Plan-mi nus-hase  differences  in  the  velocity  con¬ 
stants  for  the  four  steady-state  tests  and  summary  statistics  arc  given  in 
Tables  10A-10D. 

64.  The  inertial  and  frictional  forces  inherent  in  the  velocity  meters 
introduce  relatively  large  instrument  errors  at  lew  current  velocities; 
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however ,  as  the  cunent  speed  increases  these  instrument  errors  decrease. 

Thus  a  higher  degree  of  confidence  can  be  placed  in  the  discrete  maximum  flood 
and  maximum  ebb  velocities  measured  in  the  model  such  that  plan-to-base  com¬ 
parisons  are  meaningful.  Each  set  of  25  hourly  velocities  contained  two  maxi¬ 
mum  flood  and  two  maximum  ebb  velocity  measurements.  These  two  maximums  were 
averaged  and  are  listed  in  Tables  11A-11I)  under  the  "MODEL  DATA"  heading. 

Also,  maximum  flood  and  ebb  velocities  computed  from  the  harmonic  Fourier 
analysis  are  given  in  Tables  11A-11D  under  the  "HARMONIC  ANALYSIS"  heading. 
Plan-minus-base  maximum  flood  and  ebb  differences  and  summary  statistics  are 
given  in  Tables  12A-12F.  Plates  1 13-144  show  the  maximum  flood  and  ebb  ve¬ 
locity  profiles  for  both  the  base  and  plan  tests. 

65.  The  concept  of  flow  predominance  is  useful  in  analyzing  the  effects 
of  density- induced  currents  on  velocities.  Consider  the  conventional  velocity 
versus  time  plot  in  Figure  18.  The  shaded  area  under  the  ebb  portion  of  the 
curve  is  divided  by  the  sum  of  the  areas  under  the  ebb  curve  and  flood  curve. 
The  result  is  called  the  flow  predominance  and  is  defined  as  the  percent  of 
the  total  flow  per  tide  cycle  that  is  moving  in  the  ebb  direction  at  a  given 
velocity  sampling  depth.  In  a  partly  mixed  estuary  such  as  Chesapeake  Bay, 
the  bottom  flow  predominance  will  be  mainly  in  the  flood  direction  within  the 
salinity  wedge,  and  the  surface  flow  predominance  will  be  in  the  ebb 


Figure  18.  Definition  sketch-  flow  predominance 
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direction.  lhe  t  low  predominance  percentages  lor  the  tour  steady-state  tests 
along  with  plan-minus-base  differences  and  summary  statistics  are  listed  in 
Tables  13A-13B. 


Velocity  Test  Results 

P 1  an- tt)-base  velocity  phase  comparisons 

6b.  Current  velocity  phase  is  the  time  (in  degrees)  ot  maximum  flood 
velocity  referenced  to  zero-lunar-hour.  A  complete  tidal  cycle  in  the  model 
is  12.4206  model  hours  (360  deg)  in  length.  The  confidence  of  velocity 
phasing  is  approximately  120  deg.  Frequency-of -occurrence  statistics  of  plan- 
minus-base  phase  differences  are  provided  in  Table  10B.  Nearly  TO  percent  of 
the  plan  and  base  values  were  within  5  deg  of  one  another,  67  percent  were 
within  10  deg,  nearly  87  percent  were  within  13  deg,  and  nearly  97  percent 
were  within  the  20-deg  onfidence  band.  Table  1 0A  shows  that  during  Test  1 
the  plan  phase  occurs  64  deg  prior  to  the  base  phase  at  sta  JG0103  at  depth 
81  ft.  A  velocity  meter  problem  was  evident  during  the  base  test  and  this 
phase  value  should  not  be  considered.  Also,  at  sta  FI10901  and  depth  33  ft 
during  Test  2  and  .it  sta  F.F.0301  and  depth  21  ft  during  Test  4  the  velocities 
were  so  low  that  t h resho 1 d-o 1 -mot l on  problems  were  experienced  resulting  in 
anomalous  phasing  Only  10  (2.0  percent)  phase  discrepancies  fall  outside  the 
20-deg  confidence  band  with  all  but  one  occurring  in  the  Klizabeth  River  at 
sta  F.H0202,  EH0201,  ind  EH0301.  It  appears  that  the  deepened  Elizabeth  River 
channel  with  its  larger  cross-sectional  area  reduced  velocities  sufficiently 
to  cause  the  velocity  meters  to  occasionally  operate  within  their  threshold- 
of-motion  ranges  resulting  in  the-  anomalous  data.  Thus  any  phase  differences 
outside  the  tZ 0  deg  confidence  band  should  be  viewed  with  reservation. 

67.  The  conclusion  is  that  channel  deepening  has  caused  only  minimal 
phasing  differences  between  base  and  plan  tests.  The  only  sampling  depth  at 
which  a  phase  difference  trend  seemed  evident  was  at  sta  LH0501,  depth  39  ft, 
which  is  at  the  confluence  of  the  Eastern  and  Southern  Branches  of  the  Eliza¬ 
beth  River.  Here  the  plan  phase  consistently  arrived  earlier  than  the  base 
phase . 

PI  an- to-base  velocity 
amplitude  comparisons 

of! .  Cuiifiit  vrlwi  ity  amplitudes  are  defined  as  one-h3lf  of  the  range 
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between  till*  max  I  iiiiuii  i  looii  ami  maximum  ebb  vt-im  ilies  i  i  urn  I'm-  losiuc 
All  23  hourly  vclm  itv  obse  rva  t  i  mis  a  re  taken  11  ito  arrmuit  in  the  amplitude 
c omput  it  ions,  so  greater  ronl  iitcine  can  he  placed  in  the  amplitude  values  than 
in  a  discrete  velocity  observation.  Frequent v-of -occurrence  statistics  of 
p  Ian-1111  nus-hase  amplitude  differences  are  given  in  Table  IOC.  fiver  27  percent 
of  the  base  and  plan  velocity  amplitudes  were  within  0.10  tps  of  one  another, 
over  hi  percent  were  within  0.20  fps,  over  75  percent  were  within  0 .  50  tps, 
and  about  84  percent  were  within  0.40  tps.  There  was  an  overall  trend  toward 
decreased  amplitudes  in  the  plan  test.  Of  the  547  amplitude  difference  values 
9 7  f28  percent)  were  greater  than  0.00  fps  while  230  (72  percent)  were  less 
than  0.00  fps.  The  overall  mean  of  the  amplitude  differences  indicated  a 
0.1 3- fps  decrease  in  the  plan  test.  This  decrease  in  velocity  amplitude  in 
the  plan  test  meets  the  expectancy  associated  with  deepened  channels  and  the 
resulting  larger  cross-sect lona 1  areas. 

69.  Anomalous  amplitude  values  are  evident  at  sta  CB0008  and  sta  JG0105 
depths  6f>  and  85  ft,  respectively,  during  base  Test  1  and  should  not  be  con¬ 
sidered  due  to  meter  instrument  errors.  The  largest  plan-to-base  decreases  in 
amplitude  are  between  0.80  and  1.00  fps  and  occur  during  Test  1  (200,000-cfs 
discharge  with  spring  tides)  at  sta  100102,  F.110202,  and  EH0203.  The  largest 
plan-to-base  increases  in  amplitude  are  0.36  fps  at  sta  CL50006,  depth  4  ft, 
during  Test  1  and  0.31  fps  at  sta  YS0001,  depth  23  ft,  during  Test  2. 

70.  By  inspection  of  Tables  1 0A ,  10B,  and  100,  trends  can  be  seen  in 
the  amplitude  changes  due  to  channel  deepening  without  regard  to  a  specific 
freshwater  discharge  or  tide  range.  On  the  average,  the  stations  on  range 
CBU0  show  plan-to-base  decreases  in  amplitude  with  the  exception  of  0150004  and 
CB0006  which  show  slight  increases.  An  overall  decrease  in  amplitude  in  the 
plan  test  is  also  evident  across  range  CB01  with  the  exception  of  the  surface 
depth  at  OB0109.  Sta  A00002  shows  an  amplitude  decrease  at  the  surface  and 
middle  while  t  fie  bottom  amplitudes  remained  essentially  the  same  in  the  plan 
tests.  Sta  YS000 1  shows  an  amplitude  decrease  at  the  surface  and  averages 
essentially  no  change  at  the  middle  and  bottom.  Sta  TS000 5  indicates  a  de¬ 
crease  in  surface  and  middle  amplitudes  and  es'-lialtv  no  change  .it  the  bot¬ 
tom.  Sta  TSOOOr>  indicates  a  decrease  at  the  surface,  no  change  in  the  middle, 
and  an  increase  at  the  bottom.  Interestingly,  the  entrance  to  the  James  Rivet 
at  range  JG01  shows  plan  decreases  in  amplitude  at  the  out -ol -channel  stations 
while  the  deep  natural  channel  sta  JG010.5  shows  an  increase  in  amplitude  it 
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.i  redistribution  of  the  flow  pattern  at  t  fie  muitli  of  t  tie  lames.  1  lie  other  Un 
lames  ranges,  JN02  anti  JG03  ,  show  s  1  l  glit  ilee  teases  in  amplilmie  essentially 
across  botft  ranges.  The  stations  in  tlie  Kli/aheth  Kiver  anil  its  tributaries 
all  show  ilei  reases  in  plan  amp  1  1 1  titles  except  the  hot  lorn  depths  at  sla  EH0501 
ami  EH0901  which  indicate  essentially  no  change  f i um  base  to  plan. 

Plan-to-base  mean  velocity  comparisons 

71.  The  mean  current  velocity  is  simple  the  mean  oi  the  25  hourly 
velocity  measurements  at  a  given  sampling  depth.  The  mean  velocity  gives  an 
indication  of  flow  predominance,  positive  lie  i  ng  flood  pt  edoin  i  nance  and  negative 
being  ebb  predominance.  The  f requeue y-of-nci urrem e  statistics  of  plan-minus- 
hase  mean  velocity  differences  are  listed  in  Table  1 01 J .  More  than  60  percent 
of  the  base  and  plan  mean  velocities  were  within  0.10  fps  of  one  another,  over 
86  percent  were  within  0.20  fps,  nearly  9r>  percent  were  w  i  t  h  i  n  0.  10  fps,  and 
over  97  percent  were  within  0.40  fps.  The  mean  of  these  difference's  for  each 
of  the  four  s t eady- s t a t e  tests  was  close  to  0.00  fps  is  was  the  overal I  mean 

of  the  entire  data  set.  This  indicates  that  channel  deepening  had  little  im¬ 
pact  on  the  mean  current  velocities  when  t fit'  entire  data  set  is  considered. 

72.  The  average  overall  p 1 an-m i  mis  -  base  mean  velocity  differences  at 
each  sampling  depth  for  the  four  tests  are  given  in  Table  101).  Inspection  of 
Table  10A  shows  that  sampling  points  at  which  the  mean  velocity  increased 
with  channel  deepening  during  all  four  test:  are: 

CB0001  Depth  4  ft 
TS0003  Depths  4  ami  4b  if 
JG0103  Depths  22  and  66  ft 
JN0204  Depth  48  ft 
EH0701  Depth  4  ft 

Sampling  points  showing  a  decrease  in  plan  mean  velocity  for  all  four  tests 
a  re : 


CB0004  Depth  4  ft 
TS0005  Depth  50  ft 
EH0202  Depth  46  ft 

All  other  stations  show  either  essentially  no  change  from  base-to-plan  condi 
tions  or  an  increase  during  one  test  and  a  decrease  in  another. 
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P lan- to-base  maximum  flood 
and  ebb  velocity  comparisons 

73.  Since  the  maximum  flood  and  ebb  velocities  are  normally  well  out¬ 

side  the  thresho Id-of -mot  ion  range,  a  high  degree  of  confidence  can  be  placed 
on  these  measurements  making  plan-to-base  comparisons  meaningful.  The  plan- 
minus-base  maximum  flood  and  ebb  velocity  differences  are  presented  in  two 
forms:  (a)  the  differences  using  discrete  model  measurements  are  given  in 

Table  12A,  and  (b)  the  differences  using  the  peaks  of  the  harmonic  cosine-fit 
curve  are  given  in  Table  12B.  Anomalous  data  at  sta  CB008,  all  depths,  and 
JG0103,  depth  83  ft,  during  Test  1  were  not  considered  in  the  analyses. 

74.  From  the  f requency-of -occurrence  statistics  in  Tables  12C-12F,  the 
following  summary  of  the  overall  data  set  shows  the  percent  occurrence  within 
certain  intervals: 

Maximum  Flood  Maximum  Ebb 

percent  percent 


Cos i ne- 

Cos i ne- 

Mode  1 

Fit 

Mode  1 

Fit 

Interval 

Data 

Data 

Dai 

L  a 

Dat  a 

Plan 

and 

base 

within 

0. 

,  10 

fps 

27.9 

29.9 

21 

.  3 

22.3 

Plan 

and 
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0. 

20 

fps 

51.5 

55.2 

42 

.9 

4b .  3 

Plan 

and 
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within 

0. 

,30 

fps 

69.3 

72.7 

61 

.7 

62.3 

Plan 

and 

base 

within 

0. 

,40 

fps 

81.9 

84.4 

73 

.8 

76.6 

When  considering  overall  averages  from  the  four  tests,  the  prevailing  trend 
indicates  that  the  maximum  flood  and  ebb  velocities  decrease  by  approx i ma t o 1 y 
0.10  fps  after  channel  deepening.  The  largest  average  change  was  during 
Test  4  (70,000-cfs  discharge  with  neap  tide  range)  where  plan  maximum  flood 
and  ebb  velocities  decreased  about  0.20  fps.  The  overall  maximum  flood  and 
ebb  velocity  magnitudes  were  reduced  slightly  because  of  the  increase  in  cross 
sectional  area  associated  with  the  deeper  plan  channels. 

75.  Plates  113-144  show  two  sets  of  maximum  flood  and  ebb  velocity  pro¬ 
files  for  each  steady-state  test  at  each  station.  One  set  is  labeled  "model 
data"  and  represents  the  actual  maximum  flood  and  ebb  velocities  measured  in 
the  model.  The  other  set  is  labeled  "harmonic  data"  and  represents  the  maxi¬ 
mum  flood  and  ebb  velocities  taken  from  the  harmonic  cosine-fit  curve  result¬ 
ing  from  the  harmonic  analysis  of  the  25  hourly  velocity  observations.  Most 
often  the  model  and  harmonic  profiles  are  in  agreement.  However,  there  are 
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times  when  the  harmonic  profile  shows  better  p 1 an-to-base  agreement  than  does 
the  model  data  profile  (e.g.,  see  the  ebb  profile  for  sta  TS0005  during  Test  2 
in  Plate  127).  Likewise,  there  are  occasions  when  the  model  data  profile 
shows  slightly  better  plan-to-base  agreement  than  does  the  harmonic  profile 
(e.g.,  see  the  ebb  profile  for  sta  CB0008  during  Test  2  in  Plate  117).  Both 
sets  of  velocity  profiles  were  given  equal  weight  and  were  used  as  references 
to  formulate  the  following  specific  comments: 


a.  The  maximum  flood  and  ebb  velocities  at  range  CBOO  for  both 
base  and  plan  conditions  were  essentially  the  same  with  the  ex¬ 
ception  of  sta  CB0002  which  shows  lower  ebb  velocities  during 
the  plan  Tests  1  and  2  (200,000-cfs  discharge),  sta  CB0004 
which  showed  increased  flood  and  ebb  during  Test  2  and  sta 
CB0008  which  showed  reduced  flood  and  ebb  during  Test  4.  De¬ 
creases  in  plan  maximum  velocities  can  be  seen  at  sta  CB0009 . 

b.  Across  range  CBO 1 ,  the  overall  trend  was  for  maximum  flood  and 
ebb  velocities  to  decrease  slightly  along  the  western  half  of 
the  range  and  remain  essentially  the  same  along  the  eastern 
half  during  the  plan  test. 

c.  At  sta  AC0002,  the  effects  of  the  new  Atlantic  Ocean  channel 
can  be  clearly  seen  in  Plate1  125.  The  deepened  channel  with 
its  larger  cross-sectional  area  causes  a  decrease  in  plan 
velocities  at  the  surface  ami  middle  depths. 

d.  At  sta  TS0003,  the  plan  maximum  ebb  velocities  tend  to  decrease 
at  the  surface  while  the  flood  velocities  show  little  change 
except  during  Tests  1  and  4  where  there  is  a  decrease  of  about 
0.40  fps  at  the  surface  and  middle  depths.  At  TS0005  ,  the  max¬ 
imum  velocity  profiles  are  essentially  the  same  from  base-to- 
plan  conditions  except  for  t  tie  middle  depth  during  Test  4  which 
becomes  more  flood  predominant. 

e.  At  the  mouth  of  the  .James  River,  sta  JG0101  and  JG0102  both  in¬ 
dicate  a  general  decrease  in  maximum  velocities  of  about  0.20 
to  1.00  fps.  Sta  JG0103  (Plate  131)  contains  anomalous  data  at 
depth  83  ft  during  base  Test  1;  however,  there  is  an  indication 
of  increased  velocities  at  bottom  depths  during  the  spring  tide 
range  (Tests  I  and  3)  while  the  entire  velocity  profile  remains 
unchanged  during  neap  tides  (Tests  2  and  4). 

f .  At  JN0204,  which  is  stationed  within  the  deepened  Norfolk 
Harbor  Channel,  a  decrease  of  about  0.40  fps  occurs  generally 
throughout  the  entire  velocity  profile  during  Test  3  while  an 
increase  of  about  0.50  fps  in  surface  ebb  is  apparent  during 
Test  1  with  a  lower  increase  for  Test  2. 

g.  Sta  JG0302  is  interesting  because  a  0.20-  to  1.00-fps  decrease 
is  evident  throughout  the  plan  flood  and  ebb  velocity  profiles 
during  Tests  3  and  4  (70,000-cfs  discharge).  However,  there  is 
an  increase  of  about  0.35  fps  throughout  the  entire  flood  and 
ebb  profile  during  Test  1.  Test  2  shows  essentially  no  change. 
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ll .  Sla  EH0202  is  another  good  example  showing  lioW  channel  deepen¬ 
ing  causes  a  decrease  in  velocities  within  the  channel  (the 
base  velocities  during  Test  2,  however,  appear  to  be  within  the 
threshold-of-motion  range  of  the  velocity  meter  and  should  not 
be  considered).  All  velocity  stations  within  the  deepened 
Elizabeth  River  channel  show  a  general  decrease  in  maximum 
flood  and  ebb  velocities  of  0.10  to  0.40  fps  throughout  the 
vertical  profile. 

Plan-to-base  flow 
predominance _ compa  r i sons 

76.  Flow  predominance  values  are  given  in  percent  of  flow  in  the  ebb 
direction  in  Table  13A.  Frequency-of-occur rence  statistics  for  plan-minus- 
base  flow  predominance  differences  are  given  in  Table  13B.  Overall,  about 
66  percent  of  the  differences  were  within  5  percent,  over  83  percent  were 
within  10  percent,  about  92  percent  were  within  15  percent,  and  nearly  97  per¬ 
cent  were  within  20  percent.  The  overall  mean  difference  was  -0.57  percent. 
Tests  1  and  2  with  the  high  200,000-cfs  discharge  show  slightly  overall  in¬ 
creased  ebb  predominance  during  the  plan  conditions  while  Tests  3  and  4  with 
the  lower  70,000-cfs  discharge  indicate  that  channel  deepening  will  decrease 
ebb  predominance  slightly.  It  should  be  mentioned,  however,  that  the  sampling 
stations  represent  a  finite  number  of  points  and  do  not  portray  the  entire 
cross  section.  Other  stations  across  the  section,  if  they  were  sampled,  might 
have  refuted  these  observations. 

77.  A  more  general  conception  of  flow  conditions  in  the  James  River  and 
Elizabeth  River  is  shown  in  Figures  19  and  20  where  surface  and  bottom  flow 
predominance  profiles  have  been  plotted  along  the  estuary  channels.  Plan- 
versus-base  changes  are  subtle  in  the  James  River,  normally  not  varying  by 
more  than  5  or  10  percent;  however,  consistent  trends  are  evident.  Plan- 
versus-base  variations  in  the  flow  predominance  profile  .ire  more  apparent  in 
the  Elizabeth  River.  However,  it  must  be  remembered  that  velocity  meters  were 
operating  near  the  threshold-of-motion  range,  so  less  confidence  can  be  placed 
in  the  Elizabeth  River  measurements.  Again,  consistent  trends  are  evident. 

78.  Saltwater  intrusion  is  important  to  estuary  sedimentation  because 
saline  water  causes  flocculation  of  suspended  clay  particles  and  density  cur¬ 
rents  tend  to  move  sediments  upstream  along  the  bottom  Thus  sediments  enter¬ 
ing  the  estuary  may  become  trapped  instead  of  moving  back  out  to  sea.  Fre¬ 
quently,  the  heaviest  shoaling  occurs  between  the  high-water  and  low-water 
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JAMES  RIVER  CHANNEL 


1 000- FT  CHANNEL  STATIONS 


Figure  20.  Flow  predominance  along  the  Elizabeth  River  channel 


positions  of  the  upstream  limit  of  salinity  intrusion.  Although  shoaling  de¬ 
pends  largely  on  the  ability  of  bottom  currents  to  move  bed  materials  upstream 
the  most  likely  region  in  thich  the  heaviest  shoaling  will  occur  is  the  reach 
bracketing  the  50  percent  value  (or  null  point)  on  the  bottom  flow  predominant 
profile  (Schultz  and  Simmons  1957).  In  both  the  Janies  and  Klizabeth  Rivers  , 
the  location  ot  this  null  point  remains  basically  the  same  from  base-to-plan 
conditions.  Hence,  channel  deepening  should  not  change  the  present  locations 
of  the  heaviest  shoaling  in  the  James  and  Elizabeth  Kivers.  On  the  other  hand 
there  are  several  isolated  instances  where  the  plan  condition  reversed  the  di¬ 
rection  of  bottom  flow  predominance.  This  should  indicate  a  significant 
change  in  local  shoaling  during  those  particular  flow  and  tide  conditions. 
However,  this  model  study  cannot  directly  predict  whether  the  magnitude  of 
shoaling  will  change  with  channel  deepening.  It  is  important  to  understand 
that  the  locations  of  these  heavy  shoaling  areas  cannot  be  pinpointed  exactly 
in  the  model  and  are  only  "best  guess"  interpretations  ot  model  velocity  mea¬ 
surements.  In  order  to  determine  shoaling  locations  more  accurately,  a  more 
dense  velocity  sampling  network  would  be  necessary.  The  velocity  sampling 
scheme  used  in  this  study  was  primarily  determined  by  the  re»pi  i rement s  of  the 
numerical  hydrodynamic  models  that  will  perform  a  much  more  rigorous  sedimenta 
tion  and  shoaling  analysis  than  was  presented  herein. 

Tide  and  Velocity  Summary 

79.  Any  changes  in  tide  elevations,  amplitudes,  and  phasing  due  to  chat; 
net  deepening  were  sufficiently  small  that  they  remained  undetectable  with  the 
measurement  techniques  used  at  the  hydraulic  model. 

80.  Several  subtle  velocity  variations  in  the  model  tests  were  ap¬ 
parently  due  to  channel  deepening.  An  overall  average  decrease  in  velocity 
amplitude  ol  about  0.1J  fps  was  observed  in  the  plan  test .  This  decrease  is 
consistent  with  the  principle  of  continuity  which  states  that  for  a  given  flow 
rate  the  velocity  will  vary  inversely  with  the  cross-sectional  area.  Simi¬ 
larly,  the  maximum  flood  and  ebb  velocities  tend  to  decrease  by  approximately 
0.10  fps  on  the  average  during  the  plan  test.  During  the  high  discharge  tests 
the  overall  ebb  predominance  increased  slightly.  During  the  70,000-cts  dis¬ 
charge  tests,  which  represent  the  long-term  average  annual  t 1 ow  into  Chesa¬ 
peake  Bay,  the  overall  ebb  predominance  decreased  slightly.  But  it  must  be 
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were  located  in  the  channels,  so  they  would  teel  the  impart  of  salinity  intru¬ 
sion  more  than  the  sha 1 1 o  (-water  stations  and  would  also  tend  to  bias  the  over¬ 
all  flow  predominance  mean  differences  referred  to  m  Table  1JB.  Hence,  no 
firm  conclusions  concerning  salinity  intrusion  into  the  estuary  can  be  made 
from  the  velocity  data.  Deepening  the  Norfolk  channels  should  not  affect  the 
present  location  ot  heaviest  shoaling  in  the  James  and  Elizabeth  channels, 
based  on  the  limited  number  of  sampling  stations  used  in  this  study. 


47 


PART  IV:  DYNAMIC  SAM NI TV  TKSTING 

Test  Conditions 

81.  Dynamic  salinity  testing  for  the  Norfolk  Harbor  study  was  designed 
so  that  model  response  to  base  and  plan  geometries  could  be  simulated  for  a 
range  ot  naturally  occurring  prototype  conditions.  The  salinity  structure  of 
the  prototype,  particularly  in  the  lower  bay  and  Janies  River  area,  is  sensi¬ 
tive  to  a  number  of  boundary  conditions  that  can  be  simulated  in  the  model . 

The  two  most  important  conditions  are  the  highly  variable  nature  of  treshwute. 
input  and  the  cyclical  neap-t o-spr ing  variations  in  the  tides. 

82.  The  effect  of  freshwater  input  on  partially  mixed  estuaries  is  well 
documented  in  the  literature.  Simply  stated,  it  has  been  shown  in  numerous 
estuaries  that  increasing  freshwater  input  upstream  tends  to  stratify  the 
estuary  given  constant  tidal  conditions,  while  decreasing  inflows  provides  i 
homogenizing  effect.  It  the  degree  ot  stratification  is  not  altered  drasti¬ 
cally,  an  increase  ill  freshwater  inflow  reduces  the  upstream  extent  ot  salt¬ 
water  intrusion  and  vice'  versa.  Tidal  influence  on  salinity  structure  lias 
also  been  studied  m  the  past  but  generally  in  a  steady-stale  tashiou  by  ob¬ 
serving  different  constant  amplitude  tides  that  typify  maximum  and  minimum 
variations  in  the  tidal  record. 

83.  However,  the  prototype  tide  record  in  ('hesapeake  Bay,  as  well  as  it. 
other  estuaries,  does  not  consist  of  a  series  of  constant  amplitude  tides. 
Tides  in  the  Chesapeake  show  strong  neap  and  spring  variations  as  well  as 
smaller  semidiurnal  differences  in  tide  heights.  The  effect  ot  this  tidi! 
variability  on  salinity  structure  was  thought  to  be  important  tor  some  time 
but  there  has  been  very  little  substantiation  in  the  literature.  The'  most 
notable  prototype  studies  to  date  discussing  this  phenomenon  are  1 rom  Haas 
(1977)  and  Allen  et  al.  (1980). 

84.  Haas'  study  area  was  the  lower  York  and  Rappahannock  Rivers  in 
Virginia  where  he  found  strong  correlations  between  salinity  stratification 
and  the  neap-spring  tide  cycle.  Allen's  studies  were  primarily  confined  to 
macrotidal  estuaries  in  France,  but  his  observations  do  support  the'  importance 
of  the  neap-spring  tide  cycle  on  salinity  structure. 

83.  There  have  been  model  studies  that  have  confirmed  the  importance  of 
the  neap-spring  cycle  on  salinity  Physical  model  studies  by  Richards  and 


48 


Giilbi\iiulsf  u  (.1982)  a  tub  bitatuit  and  Gulbiandsen  fl982)  doc  1 1 1 Ti  tn  t  L  he  cxislcru  '  ’  i 
structural  salinity  changes  at  various  stations  throughout  Chesapeake  Bay. 
These  changes  seem  to  he  strongest  in  the  lower  hay  and  James  River  areas. 


8b.  Ample  evidence  exists,  therefore,  to  suggest  that  dynamic  tidal  and 
freshwater  inputs  are  necessary  to  approximate  prototype  responses  to  base  and 
plan  geometries  in  the  Norfolk  Harbor  areas.  The  Norfolk  Harbor  study  simu¬ 
lated  variable  prototype  conditions  by  using  a  2-1/2-vear  weekly  stepped  hydro 
graph  and  a  repetitive  28-lunar-day  (lib  cycle)  variable  tide.  The  source 
salinity  for  the  model  study  for  ease  ot  duplication  in  both  tests  was  held 
constant  at  32.8  ppt  which  is  appropriate  for  the  historical  drought  hydro¬ 
graph  used.  A  more  detailed  discussion  ot  test  boundary  conditions  follows. 

T  l  des 

87.  Source  (.ocean)  tides  used  in  the  Norfolk  Harbor  study  consisted  of 
a  repetitive,  28- lunar-dav ,  lib-cycle  tide  sequence.  The  12  harmonic  con¬ 
stituents  used  to  construct  the  tide  were  based  on  coefficients  used  by  NuAA 
to  predict  tides  at  Old  Point  Comfort,  Virginia.  Once  the  desired  tide  for 
Old  Point  Comfort  had  been  determined,  it  was  necessary  to  adjust  the  ampli¬ 
tudes  and  phases  to  obtain  values  for  the  ocean  source  tide  generated  some 
30  miles  offshore.  'Hie  ocean  tide  is  based  on  the  equation: 
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A  =  mean  water  level  =  -0.20  ft  N0.V0 
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a  -  constituent  amplitude,  ft 
t  =  t ime ,  h  r 

F  =  constituent  period,  hr 
tt>  =  constituent  phase,  rad 
The  following  are  the  constituent  values: 
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(example  tide  ranges  for  the  source  tide  include  4.8  ft  for  tide  1,  the  larges 
spring  tide,  and  2.6  ft  for  tide  48,  a  smaller  neap  tide. 

88.  Figure  21  is  a  graphic  representation  of  the  56-cvc ie  tide  compute 
from  the  above  constituents.  It  contains  the  spring  and  neap  tide  variations 
and  the  semidiurnal  inequalities  noted  in  prototype  tide  records.  As  in  the 
prototype,  the  tide  is  propagated  up  the  hay  and  reaches  Reedy  Point,  Delawor 
at  the  Delaware  River  end  of  the*  t.6D  Cana!  approx i ma t e  1  y  16.16  hr  iatet  . 

84.  Prior  to  the  start  <>!  the  test,  a  hriet  lead-in  period  of  constant, 
tide  amplitude  is  necessary  to  adjust  the  ocean  tide  control  amplifier.  The 
range  of  tins  lead-in  tide  reflects  t he  values  of  tide  1.  Once  the  tide  amp  1 
fiet  had  been  calibrated  for  the  lead-in  tide,  the  variable  tide  was  started 
and  along  with  a  constant  discharge  condition  the'  model  was  hi  ought  to  a  stat 
i  stable  density  equilibrium.  l’pon  reaching  density  equilibrium,  the  test 
hydrograph  was  started  and  variable  tides  were  maintained  through vat  the  test 

!()  .  1  he  strict  control  of  boundarv  and  initial  conditions  it.  a  phys  i  Ca 

model  test  is  of  prim.n  y  importance  if  any  degree  of  confidence  is  to  be  at¬ 
tached  to  the  results.  [his  is  especially  true  for  base  versus  plan  testing 
where  model  response  should  be  related  entirely  to  designed  base  versus  plan 
changes  fe.g. ,  deepened  channels)  rather  than  any  fluctuations  m  boundary 
cond i t ions  . 

9  1  .  To  ensure  source*  tide  integrity  in  both  base*  and  plan  tests,  the 
ocean  tide  is  monitored  continuously  by  a  WLD  that  is  hard-wired  to  the  com¬ 
puter  and  a  strip  .'hart  recorder  for  real-time  observation.  In  addition, 
manual  point  gage  measurements  were  taken  twice  per  56-cycle  period  on  tides 
and  8  and  on  tides  35  and  36  to  confirm  water  levels  indicated  by  electronic 
measu  remen t  s . 

92.  The  Delaware*  Bay  source  tide;  was  not  used  in  these  tests  tor 
two  reasons.  First,  available  prototype  data  are  inadequate  to  define  the 
amplitudes  and  periods  of  the  source  tide  and  sal  ini  tv  under  variable  tidal 
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ure  21.  Ore  .in  smirm*  variable  tide  show!  n>»,  tides  1,  10,  ,?H,  ant 
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model  had  shown  ih.it  the  hydrodynamics  ot  the  ChD  Lanai  are  very  sensitive, 
particularly  to  variations  in  mean  water-surface  elevation,  so  th.it  even  minor 
discrepancies  m  boundary  control  of  water-surface  elevations  have  significant 
impact  on  canal  hydrodynamics  and  thus  on  salinities  in  tile  Upper  Bay.  Sinie 
boundary  control  for  the  source  tide  in  Delaware  Bay  was  no'  capable  of  pr  - 
venting  small  discrepancies  in  water-surface  elevations,  it  was  decided  not  to 
reproduce  the  source  tide  tor  th  se  tests  so  that  any  changes  in  Upper  Bay 
salinities  from  the  base  test  to  the  future  test  would  not  he  erroneously  af¬ 
fected  by  possible  discrepancies  in  boundary  control. 

Freshwater  inflow 

‘1.1.  The  Norfolk  Harbor  test  hydrograph  consisted  ot  two  parts.  First, 
there  was  a  lead-in  period  (weeks  I  -  1  5 )  ot  constant  discharge  that  ret  lei  ted 
an  average  total  bay  discharge  of  70,000  cfs.  The  second  period  (week?.  10- 1  it.) 
depicted  weekly  averaged  prototype  conditions  between  24  Nay  l'jof  and.  17  August 
1905.  Hydrograph  plots  for  the  total  bay  discharge  along  with  the  live  rivers 
closest  to  the  subject  area  are  shown  in  Hlates  145-150. 

94.  The  15-week  constant  discharge  period  was  used  t . .  observe  model 
stability  prior  t<.  starting  the  historical  hydrograph  ami  salinity  sampling. 

The  ensuing  historical  hydrograph  was  used  because  it  was  an  avail  ible  data 

set  and  had  been  used  largely  on  other  Chesapeake  Bay  Mode i  tests.  As  a  whole, 
it  reflects  a  period  of  relatively  low  to  drought  conditions  on  the  hay.  flu  s 
should  be  considered  when  analyzing  the  data.  The  hvdrograph  does,  however, 
nave  sufficiently  seasonal  variations  in  discharge  to  he  useful. 

Source  salinity 

95.  Source  salinity  is  defined  as  the  salinity  of  supply  w.itei  as  it 
enters  the  headhay  at  the  ocean.  This  definition  was  chosen  to  facilitate  sam¬ 
pling  and  control  of  salinity.  During  varying  freshwater  discharges,  the  mix¬ 
ing  of  fresh  and  salt  water  in  the  model  ocean  causes  discontinuous  and  vari¬ 
able  salinities  in  the  vicinity  of  the  headhay.  Sampling  in  this  area  can 
yield  large  changes  in  salinity  over  relatively  short  periods  of  time,  espe¬ 
cially  during  high  freshwater  discharge  periods.  The  supply  sump  is  a  large 
and  nearly  homogeneous  volume  of  water  (analogous  to  the  prototype  ocean)  that 
reacts  slowly  to  a  spiked  hydrograph;  thus  it  can  he  measured  less  often  and 
with  greater  confidence  in  a  single  set  of  measurements.  Salinity  samples  were 
taken  hourly  (4.17  prototype  days)  in  the  headhay,  in  the  return  sump,  and  in 


the  supply  sump  throughout  eai  h  dyuaiitii  salinity  test.  Testing  ”1  tin-  llu  tr 
areas  together  helps  project  trends  in  the  source  salinity,  am)  thus  large 
drops  ill  source  salmttv  run  he  prevented  bv  adding  brine  to  the  return  sump 
when  it  shows  a  decline.  Source  salinity  was  kept  to  within  £0 .  r>  ppt  of  t  hi¬ 
des  i  red  22.8  ppt  throughout  both  base  and  plan  tests. 

9t>  .  Sump  control  in  general  was  considered  good  for  both  tests.  There 
were  short-term  deviations  caused  hy  response  to  freshwater  inflow  fluctuations 
in  both  tests  and  did  not  create  perceptible  discrepancies  in  the  hay  salini¬ 
ties.  The  supply  sump  source  salinity  was  checked  a  total  of  26b  times  coin¬ 
cident  to  both  base  and  plan  tests  from  lunar  day  -146  to  lunar  day  ‘Hi.  The 
base  and  plan  source  salinities  were  within  ±0.2  ppt  of  each  other  78.9  percent 
of  the  time,  and  within  ±0.8  ppt  97.0  percent  ol  the  time.  Only  8  of  the  269 
comparisons  were  outside  the  ±0.8  ppt  range  and  they  all  occurred  during  the 
lead-in  period  from  lunar  day  -29  to  lunar  day  2.  The  largest  plan-minus-base 
difference  in  sump  source  salinity  was  +0.8  ppt  at  lunar  day  -9.  The  overall 
base  test  average  source  salinity  was  22.49  ppt  and  the  plan  test  average  was 
12.80  ppt. 

Sa  1  i n 1 1 y  samp  1 i ng 

97.  It  was  necessary  during  testing  to  monitor  the  distribution  of 
salinity  within  the  bay  to  ensure  that  general  expected  patterns  were  develop¬ 
ing  and  that  anomalies  in  the  hay  were  explainable  within  the  context  of  the 
test  In  order  to  accomplish  this  task,  a  set  of  19  stability  monitoring  sta¬ 
tions  throughout  Chesapeake  Bay  and  the  major  tributaries  was  selected  to  give 
i  general  picture  of  havwide  salinity  structure.  These  salinity  monitoring 
stations  were  not  a  primary  part  of  the  comprehensive  testing  program  but  were 
selected  mainly  for  comparison  with  other  studies  performed  at  the  Chesapeake 
Bay  model.  The  salinity  monitoring  stations  were  used  to  determine  when  the 
model  had  reached  the  dynamic  salinity  equilibrium  required  for  test  initia¬ 
tion  and  were  a  real-time  test  monitoring  aid.  Samples  were  taken  at  slack 
after  flood  at  each  of  the  19  salinity  monitoring  stations  on  every  tide  2  and 
every  tide  20  throughout  the  tests.  The  samples  were  analyzed  immediately  aim 
hand-plotted  to  facilitate  rapid  analysis.  No  abnormalities  in  density  struc¬ 
ture  were  observed  in  either  test  and,  consequently,  confidence  in  their  com- 
pa rah i 1 i t y  is  high. 

98.  Salinities  in  the  bay  were  sampled  at  individual  stations  at  from 
station  (Figures  22  and  22  and  Table  14)  Approximately 


one  to  five  depths  per 


Figure  22.  Salinity  station  locations  within  studv  area 
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Figure  23.  Map  showing  all  of  the  salinity  stations 


500  samples  were  taken  at  each  sampled  tide  The  sampled  tides  were  1,  10,  28, 
and  48  at  slack  after  flood  which  represent  high-spring,  mean,  low-spring,  and 
neap  tides,  respectively  (Figure  21).  In  all,  57,041  samples  were  collected  in 
the  plan  test  and  55,331  were  collected  in  the  base  test.  The  number  ot  addi¬ 
tional  samples  for  the  plan  test  was  due  to  additional  sampling  in  the  deepened 
ctianne  1  . 

Sa  1  inity  Data_Ana  ljrs  i  s 

99.  The  purpose  of  the  dynamic  salinity  tests  was  to  determine  what 
impact  the  deepened  channels  would  have  on  the  salinity  distribution  within 
the  Chesapeake  Ray  system.  Salinity  sampling  stations  were  positioned  at  193 
locations  throughout  the  model  so  that  the  effects  on  the  entire  system  could 
be  determined.  The  majority  of  the  stations  were  located  near  the  project 
area  in  the  lower  bay  and  James  River  with  care  taken  in  the  selection  process 
so  that  biologically  sensitive  areas  would  be  covered.  Locations  of  the  sensi¬ 
tive  areas  were  provided  by  the  Virginia  Institute  of  Marine  Science  (VIMS) 
along  with  suggestions  for  nearhy  sampling  locations.  Where  physically  pos¬ 
sible,  these  locations  were  used. 

100.  Performing  a  detailed  data  analysis  on  193  stations  and  presenting 
the  results  m  a  technical  report  is  ail  arduous  and  unnecessary  exercise. 

Uata  from  all  19  1  stations  we  re  inspected  for  base- to-pl an  salinity  differ¬ 
ences  and  a  mote  reasonable  number  of  stations  were  selected  tor  a  detailed 
analysis  and  presentation  in  this  report.  Data  from  65  stations  are  presentee) 
in  this  report  along  with  various  giaphical  plots  that  were  used  in  the  data 
analysis.  The  stat  tops  chosen  for  in<  fusion  in  this  report  provide  an  ade¬ 
quate  coverage  ot  the  entire  bay  w r t h  mo i e  stations  selected  from  areas  which 
showed  notable  base  -  t  o  -  p  I  an  s.ilnutv  changes. 

101.  Sa  1  i  in  t  v  da  t  a  w.-rc  p  i  i  m.i  r  i  I  v  ana  1  v/.ed  by  using  several  graph  i  ca  1 

l  ec  hn  i  cpie  s  which  lie  included  .  r  c  I  ti  i  s  report  First,  the*  data  were  displayed 
from  both  tests  in  inultide  th  t  i  me  -  li  i  s  t  o  r  i  es  beneath  a  total  hay  discharge 
freshwater  hydrogtaph  (Plates  151-3151  Ihe  time-histories  display  the  entire 
data  set  at  65  stations  and  are  useful  t.-i  observing  changes  in  the  vet t leal 
salinity  a i  in  lure  resulting  from  hydrographic  and  tidal  variations. 

102.  (outlined  in  the  same  plates  are  displays  of  p  1  an -in  i  nns -base  dif¬ 
ferences  at  surface  and  bottom  and  p  1  ati-mi  uus-base  depth-averaged  cl  i  f  f  erences  , 
hot  h  plotted  against  the  same  time  scale  When  a  data  point  was  missing  from 
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one  or  both  tests  no  difference  was  listed,  thereby  resulting  in  discontinui¬ 
ties  in  the  plotted  curves.  At  stations  within  the  deepened  channels,  an  ad¬ 
ditional  depth  was  included  in  the  plan  test  showing  the  deeper  bottom.  Bot¬ 
tom  differences  were  calculated  using  the  deepest  depth  at  a  station  for  each 
test . 

103.  Salinity  depth-averages  were  computed  to  determine  the  net  salin¬ 
ity  change  integrated  over  the  water  column  at  various  locations  in  the  bay 
resulting  from  channel  deepening.  It  should  be  firmly  stated  that  one  should 
not  view  depth-averaged  differences  as  the  single  most  important  indicator  of 
base-to-plau  salinity  differences.  This  method  does  not  satisfactorily  con¬ 
sider  increased  vertical  stratification  that  often  accompanies  channel  deepen¬ 
ing.  For  a  bottom  dwelling  organism,  this  is  a  crucial  consideration.  Base 
and  plan  depth-averages  were  computed  as  shown  in  Figure  24.  There  are  some 
errors  that  are  intrinsic  to  this  method.  Salinities  at  the  surface  and  bot¬ 
tom  boundaries  are  computed  by  extrapolation  from  the  two  adjacent  depths. 

This  is  not  accurate  in  any  rigorous  sense  but  it  is  a  reasonable  approxima¬ 
tion  that  is  applied  in  a  similar  fashion  in  both  base  anil  plan  tests.  When 
differences  are  computed,  the  effects  of  this  error  are  minimized.  In  order 
to  avoid  the  introduction  of  additional  errors,  the  depth-averaging  algorithm 
only  considered  those  stations  having  a  complete  set  of  salinities  for  both 
base  and  plan  tests  at  a  given  point  in  time.  In  other  words,  it  a  salinity 
value  was  missing  from  a  certain  depth  at  a  given  point  in  time,  no  depth- 
averaged  salinity  difference  was  computed.  This  accounts  for  the  gaps  in  the 
depth-averaged  difference  plots  (e.g.,  Plate  132). 

104.  The  previously  mentioned  methods  are  adequate  for  observing  base- 
to-plan  differences  at  one-  station  over  time  but  do  not  allow  a  picture  ot 
salinity  changes  in  a  longitudinal  or  transverse  sense.  For  this  reason, 
longitudinal  isohaline  profiles  were  drawn  using  data  sampled  1 roin  the  tour 
tides  of  the  neap-spring  cycle  during  both  high  and  low  discharge  periods. 

Sa  1  i ni_t^  Re s u  1 1 s 

Mo del  accur a cy  a nd  repe atab  i  1  i  ty 

105.  Prior  to  any  discussion  of  the  impacts  of  channel  deepening  on 
salinities  in  the  bay,  it  is  prudent  to  discuss  the  accuracy  and  repeatability 
of  salinity  data  measured  in  the  model.  The  term  accuracy  refers  to  closeness 
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Figure  24.  Definition  sketch:  depth-averaged  salinity 


of  salinity  values  measuied  in  tile  model  Lu  the  "tiue"  Values  uu  Ui  i  iti^  lit  tfu- 
prototype  given  that  both  model  ami  prototype  experience  identical  boundary 
conditions.  The  level  of  accuracy  ol  the  Chesapeake  Bay  model  will  never  be 
exactly  known  because  it  is  virtually  impossible  to  reproduce  in  the  model  all 
the  boundary  conditions  affecting  the  prototype  (wind,  rain,  ship  movement, 
Coriolis  force,  etc.).  However,  by  considering  the  major  forces  affecting  the 
hydrodynamics  of  the  bay  (variable  tide,  river  inflows,  and  ocean  source  sa¬ 
linity),  an  approximation  of  the  level  of  accuracy  is  possible. 

106.  The  model  has  undergone  two  lengthy  verification  tests.  The  ini¬ 
tial  verification  was  accomplished  in  1978  (Scheffner  et  al.  1981)  and  a  re¬ 
verification  in  1981  (Granat  et  al.  in  preparation)  following  extensive  con¬ 
crete  repairs  prior  to  the  Norfolk  Harbor  study.  Although  direct  compara¬ 
bility  between  individual  model  and  prototype  data  points  is  not  always  pos¬ 
sible  due  to  different  boundary  conditions  and  other  differences,  these  veri¬ 
fication  tests  have  shown  tint  the  model  does  a  good  job  of  reproducing  the 
seasonal  salinity  variations  resulting  from  a  dynamic  freshwater  hydrograph 
and  the  neap-spring  tide  cycle.  When  used  to  determine  base-to-pl an  salinity 
differences  under  similar  boundary  conditions,  the  model  is  accurate  to  within 
±1  ppt..  Certainly  not  every  measured  data  point  will  he  accurate  to  within 
these  limits,  hut  the  overall  accuracy  of  the  estimation  of  salinity  changes 
resulting  from  channel  deepening  will  he  on  that  order.  Further  discussion  of 
the  individual  inaccuracies  that  combine  to  form  the  ±1  ppt  accuracy  level  is 
lengthy  and  will  not  be  addressed  in  tins  report. 

107.  A  closely  related  and  equally  important  aspect  of  model  Vesting  is 
the  repeatability  (or  precision)  of  model  results.  Here  repeatability  refers 
to  the  closeness  to  one  another  of  different  sets  of  model  salinity  measure¬ 
ments  taken  under  the  same  boundary  conditions.  Ideally,  it  would  seem  ra¬ 
tional  in  both  base  and  plan  tests  to  bring  the  model  to  the  same  state  of 
dynamic  salinity  equilibrium  prior  to  initiation  of  testing  using  the  same 
start-up  procedure.  In  this  manner,  one  can  be  assured  that  any  base- to-p l an 
salinity  differences  .ire  the  result  of  channel  deepening  and  not  due  to  some 
inconsistency  in  the  model  start-up  procedure.  However,  there  is  an  interest¬ 
ing  dilemma  springing  from  this  apiproach.  Will  the  model,  having  different 
base  and  plan  geometric  configurations  with  the  deepened  channels,  be  able  to 
reach  the  same  salinity  stability  under  boundary  conditions  that  are  in  every 
other  way  equivalent?  A  close  inspection  of  base-versus -p 1  an  lead-in 


sa  l  init  it's  and  a  good  do  a  I  of  r.iodr  1  experience  are  reepii  red  in  iirdc:  to  satis¬ 
factorily  answer  that  question. 

108.  This  dilemma  is  of  particular  importance  to  the  Norfolk  salinity 
tests.  Both  the  base  and  plan  tests  were  brought  to  stability  in  a  similar 
fashion  by  initially  flooding  the  model  with  fresh  water  and  then  introducing 
saline  water  from  the  ocean  source  via  a  variable  tide  generated  at  the  ocean. 
Despite  the  desired  strictness  of  control  in  the  start-up  of  both  base  and 
plan  tests,  variations  in  the  rates  of  fresh-  and  saltwater  introduction  were* 
noticed  between  tests  which  resulted  in  different  salinity  distributions  at 
equivalent  points  in  time.  Due  to  differences  in  the  model  start-up  proce¬ 
dures  and  perhaps  to  a  much  lesser  extent  the  different  model  geometries,  the 
model  did  not  reach  salinity  stability  in  the  same  period  of  time  tor  the  plan 
test  as  in  the  base  condition.  From  past  experiences  with  the  model,  it  was 
decided  that  an  additional  day  (200  tide  cycles)  of  lead-in  time  was  needed  in 
the  plan  condition  to  obtain  a  state  ot  stability  similar  to  the  base  condi¬ 
tion.  After  the  additional  day  of  lead-in,  the  overall  plan  salinities  were 
appi ox i mat  el y  0.8  ppl  higher  than  the  |»„se  salinities.  The  primary  study  area 
in  the  lower  hay  and  lames  River  areas  had  very  similar  base  and  plan  salini¬ 
ties  at  this  time1,  while  areas  showing  the  maximum  lead-in  differences  wore 
distant  from  t  tie  study  area  in  the  middle  and  upper  hay.  By  the  end  of  test¬ 
ing,  the  overall  base- 1  o-p  1  an  salinity  difference  at.  the  stability  monitoring 
stations  was  near  0.0  ppt  .  The  not  of  feet  of  differing  salinity  beginning 
conditions  caused  by  dissimilar  base - 1 o-p 1  an  start-up  procedures  would  be  to 
slightly  overestimate  the  magnitude  ot  any  increased  salinities  in  the  plan 
test.  The  areas  showing  the  greatest  base- to-p 1  an  lead-in  differences  will  be 
discussed  in  more  detail  later  in  this  report. 

Neap-spring  salinity  variations 

109.  One  of  the  purposes  of  this  model  study  was  to  determine  what  dif¬ 
ferences  in  salinity  structure  could  he  associated  with  the  proposed  channel 
deepening  project.  Findings  from  the  study  would  ultimately  he  used  to  deter¬ 
mine  if  the  project  has  any  adverse  effects  on  the  environment,  and  if  so,  per¬ 
haps  the  benefits  of  the  project  would  he  reassessed.  Channel  deepening  often 
alters  the  salinity  structure  of  estuaries.  In  most  cases,  the  net  effect 
would  be  to  increase  salinity  intrusion  on  the  bottom  with  increased  vertical 
salinity  stratification.  Net  increases  in  salt  averaged  over  the  water  column 
are  sometimes  evident  as  well.  ft  becomes  the  task  of  local  authorities  to 


00 


ultimately  determine  whether  or  not  the  noted  changes  lould  adveiselv  affect 
the  estuary,  most  particularly  the  biota.  This  test  was  designed  to  determine 
over  tune  the  nature  and  magnitude  of  salinity  variations  caused  by  the  chan¬ 
nel  deepening  so  that  these  values  could  be  provided  to  local  authorities  for 
their  planning  purposes.  However,  belore  any  degree  of  concern  is  placed  on 
salinity  differences  reported  herein,  an  understanding  of  the  prototype 
through  its  physical  model  representation  is  needed  especially  if  there  are 
naturally  occurring  salinity  variations  which  make  the  magnitude  of  base-to- 
plan  differences  comparatively  small. 

110.  The  lower  Chesapeake  Bay  along  with  the  James,  York,  and  Rappahan¬ 
nock  Rivers  experience  tidally  induced  salinity  variations  that  have  been 
noticed  in  the  prototype  (Haas  1979)  and  confirmed  in  the  Chesapeake  Bay  model 
The  neap-spring  variations,  as  they  are  called,  occur  when  lower  amplitude 
tides  (tides  28  and  48)  allow  stratification  in  the  water  column  followed  by 
larger  range  tides  (tides  I  and  10)  which  mix  the  column  again.  Intertidal 
salinity  changes  of  as  much  as  5  to  8  ppt  are  common  at  a  single  depth  along 
with  depth-averaged  salt  variations  of  3  to  5  ppt  (Figures  25-27).  Extreme 
variations,  such  as  at  sta  JG2302  (Figure  25),  can  result  in  variations  of 
depth-averaged  and  bottom  salinity  changes  of  9  ppt  and  14  ppt,  respectively. 
Neap-spring  changes  in  the  project  area  appear  to  be  greatest  during  periods 
of  high  discharge  when  the  vertical  density  gradient  is  the  strongest.  A  deli 
nition  sketch  showing  the  relationship  between  the  neap-spring  tide  cycle  and 
salinity  variations  is  given  in  Figure  28  (from  Richards  and  Gulhrandsen  1982) 

111.  The  time-histories  in  Plates  151-215  illustrate  that  this  phenome¬ 
non  is  present  throughout  the  model  to  various  degrees  with  the  strongest  van 
ations  located  in  the  vicinity  of  the  project  area.  From  each  reach  of  the 
deepened  channels,  a  representative  station  was  chosen  and  displayed  in  Fig¬ 
ures  25-27.  Sta  AC0002  is  located  in  the  Atlantic  Ocean  Channel,  TS0004  in 
the  Thimble  Shoal  Channel,  JN0103  in  the  Norfolk  Harbor  Channel,  and  EHObOl 
upstream  in  the  Elizabeth  River.  Sta  JG2302  in  the  James  River  was  chosen  be¬ 
cause  it  is  the  most  extreme  example  ol  variations,  and  sta  CB0404  in  the  mid¬ 
dle  bay  shows  very  little  if  any  neap-spring  variation. 

112.  The  mechanism  responsible  for  causing  neap-spring  variations  is 
not  suitably  described  in  the  literature  nor  is  it  sufficiently  understood. 
From  the  observations  of  data  from  various  model  tests  including  the  Norfolk 
Harbor  study,  it  appears  to  occur  in  areas  that  contain  large  vert i'  >1  density 
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igure  2b.  Typical  neap-spring  and  seasonal  salinity  variations 
in  the  lower  Chesapeake  Bay  and  Atlantic  Ocean 
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in  the  Elizabeth  Kiver  and  mid-bay  regions 


04 


1  a  t  i  oils 


gradients  which,  of  course,  are  noticed  in  the  lower  hay  and  Virginia  rivers. 
This  thesis  seems  to  be  supported  by  the  observations  that  when  neap-spring 
variations  occur  at  a  station  throughout  a  time-history,  they  are  greatest 
during  high  discharge  conditions  when  the  water  column  is  most  stratified 
(Figure  24). 

113.  The  major  emphasis  in  this  discussion  is  that  organisms  which  live 
in  the  project  area  should  he  adapted  to  survive  neap-spring  changes  in  salin¬ 
ity  without  significant  stress  to  their  populations.  Neap-spring  salinity 
changes  have  existed  for  ages  and  will  continue  to  do  so  unless  massive  geo¬ 
metric  changes  and/or  flow  suppression  schemes  are  implemented  by  man  or  na¬ 
ture.  Nevertheless,  the  plan  condition  in  this  test  was  analyzed  for  struc¬ 
tural  changes  in  this  neap-spring  tide  cycle. 

Discharge-induced  salinity  variations 

114.  Variations  in  the  supply  of  fresh  water  to  the  estuary  are  also  a 
source  of  naturally  occurring  salinity  variations.  Periods  of  high  discharge 
during  the  winter  and  spring  seasons  generally  freshen  the  entire  water  column 
with  an  increased  vertical  salinity  s t ra t i f i ca t i on .  With  the  approach  of  the 
summer  and  early  autumn  low-flow  periods,  the  salinity  in  the  entire  water 
column  gradually  increases  and  the  overall  stratification  decreases.  Longitu¬ 
dinal  salinity  structure  changes  induced  by  freshwater  flow  variations  are 
very  large  at  times,  especially  in  the  higher  discharge  rivers  of  the  bay. 
Smaller  discharge  rivers  also  experience  structural  changes  relar  1  to  the 
hydrograph  but  to  a  lesser  extent  due  in  part  to  the  fact  that  ...  1  flows  in 
and  out  of  the  river  are  large  in  comparison  with  freshwater  discharge. 

115.  The  James  River  shows  the  largest  seasonal  variations  in  salinity 
in  the  study  area.  It  is  one  of  the  higher  discharge  rivers  on  the  bay 
(Plate  149),  and  the  range  ol  flow  magnitudes  between  high  and  low  conditions 
is  great.  The  flushing  effect  of  a  high  discharge  rate  can  be  observed  in  the 
James  River  longitudinal  isohaline  profiles  (Plates  216  and  217).  In  general 
over  all  tidal  conditions,  the  1-ppt  isohaline  is  located  approximately 
125,000  ft  (23.5  miles)  farther  downstream  during  the  high-flow  period  than 
during  the  low-flow  period.  Similarly,  the  15-ppt  isohaline  is  about  50,000  It 
(9  miles)  farther  downstream  during  the  high-flow  conditions. 

116.  The  Elizabeth  River  can  be  classified  as  a  low  discharge  river 
(Plate  150).  Flows  in  and  out  of  the  river  are  largely  tidal  in  nature  and  as 

a  result  the  Elizabeth  is  a  rather  well-mixed  estuary  with  the  20-ppt  isohaline 
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extending  upstream  to  sta  EH1001  dutiug  the  low-flow  period  (Plate  218).  The 
Elizabeth  River  does  not  exhibit  the  well-defined  "wedge"  of  salinity  that  was 
observed  moving  upstream  and  downstream  in  the  James  River  with  changes  in 
freshwater  flow  magnitude.  However,  the  estuary  clearly  becomes  much  more 
stratified  during  the  high-flow  period  (Plate  219).  The  20-ppt  isohaline  in 
the  vicinity  of  sta  EH1001  is  replaced  by  a  6-  or  8-ppt  isohaline.  During  the 
high-flow  period,  the  surface  salinities  along  the  estuary  are  8  ppt  or  more 
fresher  than  during  the  low-flow  period,  but  bottom  salinities  are  only  about 
1  ppt  fresher. 

117.  The  Atlantic  Ocean  and  Thimble  Shoal  Channels  are  located  in  open 
areas  of  the  bay  and  ocean  and  do  not  have  a  single,  most  predominant  source 
of  fresh  water.  They  do,  however,  experience  similar  longitudinal  reactions 
to  the  variable  hydrograph  (Plates  220  and  221)  observed  in  the  James  ami 
Elizabeth  Rivers.  In  the  upper  water  column  near  sta  TS0005 ,  the  17-  to  20-ppt 
isohaline  which  is  present  during  high-flow  conditions  is  replaced  by  a  25-  to 
26-ppt  isohaline  in  the  low-flow  period.  At  the  bottom,  the  31-ppt  isohaline 
remains  in  essentially  the  same  position  during  both  high-  and  low-flow  pe¬ 
riods  at  tides  1  and  10.  However,  at  tides  28  and  48  the  31-ppt  isohaline  in¬ 
trudes  about  50,000  ft  farther  upstream  in  the  Thimble  Shoal  Channel  during 
low-flow  periods.  The  strength  of  neap-spring  variations  in  this  region  is 
dominant  as  was  shown  in  Figure  26. 

Base -to-plan  s a  1 i n lty  differences 

118.  With  an  understanding  of  the  dynamic  nature  of  the  lower  bay  under 
undeepened  existing  conditions,  we  can  now  address  changes  in  salinity  struc¬ 
ture  which  may  be  expected  with  the  proposed  channel  deepening.  Changes  that 
occurred  in  the  deepened  plan  condition  were  examined  in  four  ways:  (a)  plan- 
minus-base  differences  at  each  depth  (most  importantly  the  bottom  depth), 

(b)  plan-minus-base  differences  in  depth-averages,  (c)  channel  versus  shallow- 
water  redistributions  of  salinity,  and  (d)  plan  versus  base  changes  to  the 
neap-spring  salinity  response.  Each  of  these  will  be  discussed  separately. 

It  should  again  be  stated  in  this  discussion  of  differences  that  no  signifi¬ 
cant  changes  due  to  channel  deepening  were  noticed  in  the  model  north  of 
range  CB01.  Any  changes  that  are  displayed  in  the  plates  are  the  result  of 
boundary  and  initial  condition  problems  in  model  control  and  not  channel 
deepening . 

119.  Plan-minus-base  differences  at  each  depth  were  plotted  through 
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time  in  Plates  1  fa  1  -  2  1  c» .  It  would  he  convenient  to  generalize  on  t  ho  nature  o  f 
salinity  differences,  but  it  is  difficult  because  each  individual  station  lies 
in  a  portion  of  the  estuary  where  freshwater  and  tidal  flows  have  different 
relative  strengths  as  well  as  different  local  geometries.  Each  station  should 
he  inspected  individually  for  differences  at  depth  with  the  realization  that 
there  is  noise  in  the  data  set  so  no  particular  significance  should  be  as¬ 
signed  to  any  one  point.  Also  it  should  be  realized  that  the  degree  of  model 
density  stability  increases  with  time  so  data  collected  in  the  later  portions 
ot  the  test  are  more  realistic  in  predicting  differences.  It  is  safe  to  gen¬ 
eralize  that  channel  stations  showed  increased  bottom  salinities  with  an  over¬ 
all  increase  in  vertical  salinity  stratification.  Shallow-water  areas  and 
channels  that  were  not  deepened  showed  less  structural  variations  and  less  in¬ 
creased  salinity  than  deepened- channe  i  stations.  Figures  summarizing  those 
findings  will  be  presented  later. 

120.  Depth-averages  were  computed  to  predict  what  net  imports  or 
exports  ot  salt  could  fie  expected  in  the  project  area.  Results  from  65  study 
stations  were  grouped  into  five  categories  describing  the  average  net  increase 


dec  rease 

in  depth 

-  ave rage 

s  throughout 

the  test . 

The  f o  1  1 

owing  summarizes 

observed 

changes 

Net  Increase,  ppt 

Net  Decrease 

0.31  to 

1.01  to 

1.51  to 

ppt 

0 

.00  to  0 

.50 

1  .00 

1  .  30 

2.00 

0.00  to  0.50 

AC 000 2 

JG0601 

PG0101 

EF.0101 

EH0202 

KH0 102 

C B0 00 2 

CH01 09 

JN0104 

POO 20 2 

EHO40 1 

EH0S01 

CB0008 

CB0208 

JN0201 

RG0102 

EH 09 01 

E110601 

J NO  1 0 1 

C B0 30 3 

JN0202 

RG0202 

EH i 001 

E HO  701 

JN0103 

CR0404 

JN0203 

RG0501 

JG0203 

EH 0801 

JN0105 

CB0505 

JN0502 

TS0002 

J GO 30 2 

JN0205 

CB0611 

JN0603 

WB0201 

JN0102 

1.F0  101 

CB0705 

JN0801 

WOO  1 0 1 

JN0204 

1.S000  1 

CG0101 

MF000  1 

YG0102 

J NO 303 

I.S0002 

EH0301 

NS0101 

YG0302 

NN0001 

LS0003 

J GO  1 03 

NS0201 

YG0501 

TS0004 

1.S0004 

.JG0402 

NS0301 

JG0502 

NS0401 

121  . 

As  shown 

in  the 

preceding  list,  stations 

upstream  of  range  CB01  as 

well  as  the  upper  James  River  show  little  if  any  changes  in  their  depth- 
averaged  salinity.  Channel  stations  in  the  Elizabeth  and  lower  James  Rivers 
along  with  the  upper  end  of  Thimble  Shoal  Channel  show  noticeable  average  net 
salinity  increases  from  0.5  to  1.7  ppt  .  Some  stations,  primarily  in 
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slid  I luw*wdtei  areas,  actually  showed  depth-averaged  freshening  in  the  plan 
test  although  very  slight,  always  less  than  0.5  ppt  'esher. 

122.  The  freshening  of  shallow-water  areas  in  the  Horseshoe  Shoal  area 
is  an  example  of  a  transverse  redistribution  of  salt  in  the  cross  section  due 
to  channel  deepening.  Freshening  rarely  occurred  at  any  of  the  stations  un¬ 
less  they  were  adjacent  to  a  deepened  channel  area.  Figure  29  shows  depth- 
averaged  differences  across  the  Horseshoe  Shoal  transect  ( LS000 1 -LS0004 )  with 
sta  TS0004  (Figure  22)  being  the  deepened  station.  The  1-ppt  approximate  in¬ 
crease  in  depth-averaged  salinity  at  sta  TS0004  and  the  slight  freshening  at 
sta  LS0001,  LS0002,  LS0003,  and  LS0004  illustrate  this  point.  During  the 
first  high-flow  season,  the  redistribution  is  particularly  strong.  Other 
transects  show  a  similar  response  to  varying  degrees  in  both  the  bottom  dif¬ 
ferences  and  depth-averaged  differences  plots. 

123.  Bottom  differences  were  one  of  the  most  important  results  of  the 
study.  ft  is  on  the  bottom  where  certain  biota,  particularly  shellfish,  are 
particularly  susceptible  to  adverse  salinity  changes.  Maps  summarizing  bottom 
differences  for  the  various  tide  and  hydrograph  conditions  are  given  in 
Plates  222-23b.  These  plates  show  that  channel  deepening  increases  salinity 
along  the  bottom  of  the  channels  by  up  to  4.0  ppt,  but  shallow-water  areas  and 
undeepened  channels  show  very  little  change.  In  fact,  outside  of  the  Eliza¬ 
beth  and  lower  James  channel  areas,  increases  were  rarely  as  great  as  1  ppt  for 
any  tide  or  hydrograph  condition  tested. 

124.  The  bottom  differences  discussed  thus  far  are  viewed  at  discrete 
points  in  time  during  particular  tidal  or  hydrographic  conditions.  It  is  also 
important  to  view  the  change  in  the  bottom  variations  through  time  to  detect 
any  possible  change  to  the  neap-spring  salinity  response.  Changes  to  the  neap¬ 
spring  response  are  most  important  on  the  bottom  for  previously  stated  reasons. 
They  are  also  easily  observable  on  the  bottom. 

125.  An  inspection  of  the  base  and  plan  time-histories  in  Plates  151- 
215  will  reveal  the  magnitude  of  the  cyclic  neap-spring  response  on  the  bottom 
depths.  For  example,  in  Plate  163,  sta  F.H0102  shows  a  different  neap-spring 
response  between  base  and  plan  conditions.  The  base  contains  an  approximate 
3-  to  4-ppt  variation  throughout  the  hydrograph  while  the  plan  shows  only  a 

2-  to  3-ppt  range.  This  damping  of  the  tidally  induced  neap-spring  response  on 
the  bottom  is  prevalent  in  the  deepened  channel  areas  of  the  Elizabeth  River. 
From  sta  EH0102  upstream  to  EH080I,  the  damping  is  fairly  consistent  with  base 
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Figure  29.  Depth-averaged  salinity  redistribution  across  the  Horseshoe 
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ranges  varying  between  2  and  3  ppt  and  the  plan  between  1  and  2  ppt .  On  the 
average,  a  1-ppt  damping  was  noticed  with  the  exception  of  sta  EH0301  which 
showed  an  increased  bottom  range.  This  station,  however,  exists  in  a  portion 
of  the  river  which  was  not  deepened.  Both  base  and  plan  tests  had  their  sa¬ 
linity  probes  located  at  the  same  depths.  The  observation  of  increased  bottom 
range  in  the  plan  condition  at  sta  EH0301  is  consistent  with  the  48-ft  depth 
downstream  at  sta  EH0202  and  the  44-ft  depth  upstream  at  sta  EH0401.  Both  ex¬ 
perienced  slight  range  increases  in  what  were  the  base  bottom  depths. 

126.  Since  three  channel  sections  were  dredged  in  the  lower  bay  and 
James  River  area  it  is  interesting  that  only  the  deepened  portions  of  the 
Elizabeth  River  showed  any  distinguishable  changes  in  the  neap-spring  salinity 
response.  A  possible  explanation  for  this  observation  stems  from  the  fact 
that  each  of  the  channel  sections,  Thimble  Shoal,  lower  James  River,  and  the 
Elizabeth  River,  lie  in  radically  different  portions  of  the  estuary.  The 
James  River  is  a  high  discharge  river  with  broad  transverse  dimensions  com¬ 
pared  with  its  relatively  narrow  navigation  channel  widths.  The  Thimble  Shoal 
Channel  lies  in  the  lower  bay  where  the  significant  freshwater  influences  of 
the  James  and  other  rivers  are  felt  plus  again  it  has  a  narrow  navigation  chan¬ 
nel  in  a  broad  cross  section.  Both  are  highly  dynamic  zones  of  turbulent  mix¬ 
ing  with  strong  freshwater  and  tidal  mixing  inputs  and  strong  current  veloci¬ 
ties.  The  Elizabeth  River,  by  contrast,  is  a  narrow  river  with  negligible 
freshwater  input,  weak  current  velocities,  and  its  deepening  represents  a 
larger  portion  of  its  total  conveyance  area.  Relatively  little  mixing  is 
present  with  successive  neap  and  spring  tide  cycles  having  little  influence  on 
the  bottom  current  velocities  and  thereby  these  salinities. 

127.  The  biological  ramifications  of  a  damped  neap-spring  response  are 
worthy  of  consideration.  A  variety  of  organisms  can  tolerate  increases  or  de¬ 
creases  in  salt  concentration  provided  they  are  of  limited  duration.  The 
damped  neap-spring  response  along  with  the  slight  increase  of  bottom  salini¬ 
ties  could  place  a  number  of  species  in  a  stressful  position.  Biologists 
would  have  to  determine  the  tolerance  levels  for  each  species  versus  the 
changes  predicted  in  this  report.  Normally,  shellfish  are  of  the  most  concern 
in  deepening  studies  both  because  of  their  bands  of  salinity  tolerance  and 
those  of  their  predators.  One  heartening  fact  observed  in  viewing  the  data  is 
that  where  changes  do  occur  either  by  absolute  or  in  a  time-varying  neap¬ 
spring  fashion,  they  inevitably  occur  in  nonactive  biological  areas  in  the 
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PART  V:  CONCLUSIONS 


128.  Changes  in  tidal  elevations,  amplitudes,  and  phasing  which  may  be 
due  to  the  effects  of  channel  deepening  were  sufficiently  small  that  they  were 
undetectable  with  the  measurement  techniques  used  at  the  hydraulic  model. 

Model  measurement  techniques  are  of  sufficient  accuracy  that  significant 
changes  would  have  been  noticed;  therefore  none  are  expected. 

129.  Several  subtle  velocity  variations  in  the  model  tests  were  ap¬ 
parently  due  to  channel  deepening.  An  overall  decrease  in  velocity  amplitude 
of  about  0.13  fps  was  noticed  during  the  plan  test.  This  is  consistent  with 
the  principles  of  continuity,  but  the  magnitude  of  change  is  close  to  the  ac¬ 
curacy  limitations  of  model  instrumentation.  Slight  increases  in  flood  pre¬ 
dominance  were  noticed  under  average  inflow  conditions  indicating  perhaps  that 
salinity  intrusion  may  move  upstream  in  the  study  area.  This  observation  is 
consistent  with  the  observed  increased  salinities.  Overall  changes  in  model 
velocities  could  be  attributed  to  the  effects  of  channel  deepening  but  the 
magnitude  of  the  changes  is  barely  detectable. 

130.  Variations  in  the  model  salinity  distribution  were  noticed  that 
could  be  attributed  to  channel  deepening.  The  greatest  differences  were 
noticed  in  the  deepened  channel  areas  where  increases  in  the  bottom  salinities 
varied  between  0  and  4.0  ppt.  Channel  depth-averaged  increases  varied  between 
0  and  2  ppt.  Shallow-water  areas  near  the  deepened  channels  experienced  much 
less  of  a  salinity  increase.  At  times  there  was  actually  a  slight  freshening 
of  the  shallow  waters.  Stations  elsewhere  in  the  model  showed  modest  in¬ 
creases  in  depth-averaged  salinity  but  were  normally  less  than  0.3  ppt. 

131.  The  salinity  tests  documented  the  locations  of  stations  in  the 
study  area  that  exhibit  large  (commonly  as  great  as  .3  to  8  ppt)  salinity 
changes  due  to  the  neap-spring  tide  cycle.  The  entire  study  area  experienced 
these  variations  that  are  naturally  occurring  and  not  caused  by  channel  deepen¬ 
ing.  Channel  deepening  did,  however,  cause  a  slight  change  in  the  neap-spring 
salinity  response.  It  was  characterized  by  a  damping  of  the  cyclical  salinity 
variations  in  the  bottom  depths  of  the  Elizabeth  River.  Elsewhere,  little  if 
any  neap-spring  changes  were  detectable. 

132.  The  intention  of  this  report  and  its  accompanying  data  set  is  t<: 
furnish  biologists  and  planners  with  the  information  necessary  to  assess  the 
impacts  of  the  proposed  deepening  project  on  the  environment.  Almost  all 


deepening  projects  cause  some  degree  of  environmental  change  and  this  project 
is  no  exception.  Every  effort  has  been  made  to  quantify  the  extreme  natural 
salinity  variability  in  the  project  area  resulting  from  hydrographic  and  tidal 
variations.  Before  any  value  judgments  are  made  pertaining  to  the  benefits  of 
the  dredging  versus  any  possible  environmental  concerns,  a  good  knowledge  of 
the  region's  natural  variability  should  be  obtained.  The  variations  in  salin¬ 
ity  response  caused  by  channel  deepening  are  small  compared  with  the  varia¬ 
tions  that  occur  naturally  in  <.  le  region. 
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Sewage  Treatment  Plans  (STP)  and  Other  Diversions 


I.D. 

Name 

Tributary 

A 

James  River  STP 

James  River 

B 

Boat  Harbor  STP 

James  River 

C 

Army  Base  STP 

Elizabeth 

River 

D 

Lamberts  Point 
Outfal 1 

Eli zabeth 
River 

E 

Chesapeake  - 
Elizabeth  STP 

Little 

Creek 

- 

Surry  Nuclear 
Cooling  Intake 

James  River 

- 

Surry  Nuclear 
Cooling  Outfal 1 

James  River 

Latitude 

Longi tude 

Discharge 

cfs 

37°^4' 35“ 

76°32 ' 20" 

26 

36°37‘ 30" 

76°24‘ 48" 

34 

36°S5‘20" 

76°20 ' 1 0" 

23 

36°52 1 58" 

76°24 ' 00" 

50 

36°56 ' 35" 

76° 1 0 ' 1 7" 

40 

37°09’ 33" 

76°40' 1  51' 

3740* 

37°1 0 ' 09" 

76  42' 11" 

3740* 

TOTAL  173* 


*The  discharge  of  3740  cfs  for  the  Surry  Nuclear  Power  Plant  cooling 
diversion  was  not  added  to  the  model,  but  was  instead  circulated 
from  the  downstream  side  of  Hog  Island  to  the  upstream  side  to  simu¬ 
late  prototype  conditions.  Hence,  the  3740  cfs  is  not  included  in 
the  total. 


Table  2 


Freshwater  Discharges 
Steady-State  Velocity  and  Tide  Tests 


Inflow 

Number 

Tributary 

Di  scharge 

Tests  1  &  2 

,  cfs 

Tests  3  &  4 

1 

Nansemond  River 

1  ,063 

372 

2 

Chickahominy  River 

826 

289 

3 

Appamattox  River 

2,763 

967 

4 

James  River 

20,711 

7,249 

5 

York  River 

7,597 

2,659 

6 

Rappahannock  River 

8,120 

2,842 

7 

Wicomico  River  (Potomac) 

1,177 

412 

8 

Occoquan  Creek 

6,771 

2,370 

9 

Anacostia  River 

1  ,663 

532 

10 

Potomac  River 

21  ,997 

7,699 

11 

Patuxent  River 

2.517 

881 

12 

Severn  River 

660 

231 

13 

Patapsco  River 

1  ,751 

6 1  3 

14 

Gunpowder  River 

2  ,201 

802 

5  &  22 

Susquehanna  River 

106.335 

37,217 

16 

Bohemia  River 

1,103 

386 

17 

Chester  River 

1  ,434 

502 

18 

Wye  River 

543 

190 

19 

Choptank  River 

3,334 

817 

20 

Nanticoke  River 

4,636 

1  ,619 

21 

Pocomoke  River 

2,849 

9r  7 

23 

Elizabeth  River 

869 

3  14 

Total  Bay  Discharge 


200,000 


70,  >00 


Table  3 


Tide  Sta_t  i on  Locations 
S tea idy-State  Tide  Tests 


Station 

Body  of  Water 

Latitude 

Lonqi tude 

OCEAN 

Atlantic  Ocean 

36°58‘ 

54" 

75°33‘ 27" 

Sta  3 

James  River  (Old  Pt.  Comfort) 

37°00' 

18" 

7  6°1 8 ’ 50" 

CB0001 

Chesapeake  Bay 

36°56‘ 

17" 

76°00' 23" 

CB0004 

Chesapeake  Bay 

37°00' 

14" 

75°59 1 30" 

CBC308 

Chesapeake  Bay 

37°03 1 

13" 

75°58' 37" 

CB0101 

Chesapeake  Bay 

37°06' 

00" 

76°1 4' 36" 

CB0105 

Chesapeake  Bay 

37°08 ' 

51" 

76°08' 48" 

CB0109 

Chesapeake  Bay 

37°1 1 ' 

42" 

76°02 ' 07" 

JG0302 

James  River 

37°03 1 

21" 

76°35 ' 35" 

JN0202 

James  River 

36°55 1 

42" 

76°25' 44" 

EH0203 

Elizabeth  River 

36°53 1 

09" 

76°20'13" 

EH0901 

Elizabeth  River 

36°45 ' 

05" 

7  6° 1 7 ' 46" 

EE0301 

East  Branch  Elizabeth  River 

36°50 ' 

13" 

76°1 4' 40" 

WB0201 

West  Branch  Elizabeth  River 

36°51 1 

18" 

76°21 1 1 5" 

TS0003 

Thimble  Shoal  Channel 

36°  58 ' 

24" 

76°06 ' 30" 

TS0005 

Thimble  Shoal  Channel 

37°00' 

18" 

76°  1  3‘ 54" 

YS0001 

York  Spit  Channel 

37°02 ' 43” 

76°04' 23" 

Steady  -  S  t  a  te  T  i  de_s  :  200,000  c  f  s 
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Table  5 


Plan- 

-Minus-Base  Amplitude 

Di fferences 

Steady-State  Tide 

Tests 

Station 

Test  1 

Test  2 

Test  3 

Test  4 

Four  Test 
Average 

ft 

ft 

ft 

ft 

ft 

OCEAN 

0.02 

-0.06 

-0.04 

-0.01 

-0.02 

Sta  3 

-0.03 

0.02 

-0.01 

-0.03 

-0.01 

CB0001 

-0.01 

0.00 

0.05 

0.03 

0.02 

CB0004 

0.02 

-0.04 

0.03 

-0.05 

-0.01 

CB0008 

0.02 

-0.01 

-0.08 

-0.03 

-0.02 

CBOlOl 

0.06 

0.00 

0.00 

0.07 

0.03 

CB0105 

0.01 

-0.02 

0.00 

-0.05 

-0.01 

CB0109 

0.05 

-0.06 

-0.02 

0.03 

0.00 

JG0302 

-0.02 

-0.02 

0.01 

-0.02 

-0.01 

JN0202 

0.00 

0.01 

0.00 

0.01 

0.00 

EH0203 

0.03 

0.02 

0.02 

0.02 

0.02 

EH0901 

0.00 

-0.03 

-0.05 

-0.01 

-0.02 

EE0301 

0.00 

-0.01 

-0.04 

-0.03 

-0.02 

WB0201 

0.03 

0.0(7 

-0.01 

0.02 

0.01 

TS0003 

0.02 

-0.04 

-0.06 

0.02 

-0.01 

TS0005 

0.04 

N/A 

N/A 

-0.03 

0.00 

YSOOOl 

0.06 

-0.02 

-0.06 

0.00 

0.00 

Mean 

0.018 

-0.013 

-0.015 

-0.004 

-0.003 

Std  Dev 

0.027 

0.023 

0.037 

0.033 

0.032 

Plan-Minus-Base  Mean  Water  Level  Differences 
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CB0001 
CB0002 
CB0004 
CB0006 
CB0008 
CB0009 
CB01 01 
CB0103 
CB0105 
CB0107 
CB0109 
CB0110 
AC0002 
TS0003 
TS0005 
YS0001 
uGOlOl 
JG01 02 
JG0103 
JG0311 
JG0302 
JG0321 
JN0202 
JN0203 
JN0204 
EH0202 
EH0203 
EHQ501 
EH0701 
EH090I 
EE0301 

WB0201 


Table  8 

Veloci ty  Station  Locations 
Steady-State  Velocity  Tests 


Body  of  Water _ 

Chesapeake  Bay 
Chesapeake  Bay 
Chesapeake  Bay 
Chesapeake  Bay 
Chesapeake  Bay 
Chesapeake  Bay 
Chesapeake  Bay 
Chesapeake  Bay 
Chesapeake  Bay 
Chesapeake  Bay 
Chesapeake  Bay 
Chesapeake  Bay 
Atlantic  Ocean  Channel 
Thimble  Shoal  Channel 
Thimble  Shoal  Channel 
York  Spit  Channel 
James  River 
James  River 
James  River 
James  River 
James  Ri  v.er 
James  River 
James  River 
James  River 
James  River 
Elizabeth  River 
Elizabeth  River 
Elizabeth  River 
Elizabeth  River 
Elizabeth  River 
East  Branch  Elizabeth 
Ri  ver 

West  Branch  Elizabeth 
Ri  ver 


Sampling  Depths 


(ft  below  msl ) 

Latitude 

Longitude 

4, 

18, 

32 

36°56 ' 

12" 

76°00' 

18" 

4, 

20,35 

,54,71 

36°57' 

'30" 

76°00' 

‘05" 

4, 

16, 

28 

37°00 

'06" 

75°59' 

'30" 

4, 

15 

37°02' 

'00" 

75°59' 

'06" 

4, 

24, 

43 

37°03 ’ 

'18" 

75°58‘ 

'48" 

4, 

16 

37°04' 

37" 

7  5°58 ' 

'36" 

4, 

13 

37°05' 

‘58" 

76n14' 

‘43" 

4, 

14, 

24 

37°07' 

'05" 

76°12‘ 

'15" 

4, 

20, 

37 

37°08' 

36" 

76°08‘ 

59" 

4, 

15, 

27 

37°10' 

'36" 

76°04 1 

'37" 

4, 

39, 

75 

37°11' 

'48" 

76°02 1 

'01" 

4, 

15 

37°1 2 ' 

08" 

76°01  1 

32" 

4, 

27,49 

,54* 

36°54' 

06" 

75°54 1 

33" 

4, 

26,56 

,51* 

36°58' 

‘21" 

76°06' 

‘39" 

4, 

27,50 

,55* 

37°00' 

'24" 

76°14‘ 

'39" 

4, 

25, 

45 

37°02’ 

33" 

76°04' 

'31" 

4, 

11 

36°58‘ 

'54" 

76n17' 

'36" 

4, 

26, 

48 

36°59' 

’36" 

76°17' 

'54" 

4, 

22,44 

,66,83 

37°00' 

‘00" 

76°18‘ 

'12" 

4, 

11 

37°03 1 

‘06" 

76^36 1 

'16" 

4, 

16, 

27 

37°03' 

‘30" 

76n35' 

36" 

4 

37°04' 

'll" 

76°36‘ 

'15" 

4, 

12 

36°55 ' 

’23" 

76°25 1 

46" 

4, 

20 

36°56' 

'24" 

76^25 1 

22" 

4, 

26,48 

',55* 

36°57' 

04" 

76°25‘ 

07" 

4, 

25,46 

sc* 

y  J 

36  53 1 

15" 

76°20' 

04" 

4, 

22, 

40 

36^53 1 

14" 

76^20 ' 

16" 

4, 

22,39 

,45* 

36  50  1 

34" 

76°1 7 ' 

41" 

4, 

19,35 

,40* 

36n47' 

05" 

76°18' 

13" 

4, 

19, 

35 

36°45 ' 

18" 

76°1 7 ' 

44" 

4, 

21 

36°50' 

13" 

7  6° 1 4  ’ 

48" 

4, 

13 

36°51 1 

14" 

76°21  1 

20" 

Additional  depths  samples  during  plan  test  only. 
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NUMBER  IN  SAMPLE  84  83  83  88  348 
MEAN  DIFFERENCE  (FPS)  -0.13  -0.03  -0.33  -0.17  -3.10 
STANDARD  DEVIATION  ( FPS J  0.34  0.27  0.2^  0.28  0.29 


HZ?  IN  SmI'IPLE 
n  DIFFERENCE  (FP5) 
M'fiPD  DEVIATION  ( FP3) 


NUMBER  IN  SAMPLE  6-!  88  88  88  340 
MEAN  DIFFERENCE  rcPS>  -0.13  -0.03  -0.11  -0.18  -0.11 
STANDARD  DEVIATION  (EPS)  0.31  0.28  0.24  8.27  0.28 
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Table  14 


SAL  I N I TY  STATION  LQCAT I ONS 
DYNAMIC  HYDROGRAPH  TESTS 


Sampling  Depths 


Station 

Body  of  Water 

ft  below  msl 

La ti tude 

Longi tude 

AC0001 

Atlantic  Channel 

4,14,24,50 

35°49  1  54" 

75°48 ' 52 1 

AC0002 

Atlantic  Channel 

4,14,24,50,57* 

36°54 1 06" 

75°  5  4 ' 33" 

AC0003 

Atlantic  Channel 

4,14,24,50,57* 

36°55 1 10" 

75°56 '12" 

B GO  101 

Back  River 

(Vi rqinia) 

3,10 

37°06 1 36" 

76° 17' 30" 

CB0001 

Chesapeake 

Bay 

5,13,23,33 

36°56 1 12" 

76°00' 18" 

C30002 

Chesapeake 

Bay 

4,12,22,32,68 

36°57'30" 

"6°00'05" 

CB0003 

Chesapeake 

Bay 

4,22,32 

36°58 1 48" 

73°59 ’ 43" 

CB0004 

Chesapeake 

Bay 

4,22,30 

37°00 1 06" 

7559 ' 30" 

CB0005 

Chesapeake 

Bay 

4,12,17 

37°01 '12" 

75° 59  IS" 

CB0006 

Chesapeake 

Bay 

5,12 

37°02'00" 

75c  59 1 06 " 

CB0007 

Chesapeake 

Bay 

5,13,18 

37°02 1 54" 

7  5°5  3 ' 5  4 ' 

CB0008 

Chesapeake 

Bay 

5,23,43 

37°03' IS" 

7  5  °  5  S ‘ 43" 

CB0009 

Chesapeake 

Bay 

4,16 

37°04’ 37" 

7  5  0  5  8 ' 36 

CB0101 

Chesapeake 

Bay 

5,15 

37°05'j3" 

7  6° 1 4 ' 43" 

CB  0 102 

Chesapeake 

Bay 

4,12,13 

3  7°OC '  3 7 " 

76°’ 3'?  ' 

CB0103 

Chesapeake 

Bay 

5,15,28 

37°07 ' 05" 

76°12 ’ 15" 

CB0104 

Chesapeake 

Bay 

4,13,30 

37°07 1 37" 

76°1 1 1 07" 

CB0105 

Chesapeake 

Bay 

4,23,37 

37°08' 36" 

76  08 ' 59 

CB0106 

Chesapeake 

Bay 

4,12,22 

37°09 ' 44" 

76° 06 ' 75" 

CB0107 

Chesapeake 

Bay 

4,12,27 

37°10 ' 36" 

7CC0  *'  37" 

CB0108 

Chesapeake 

Bay 

4,14,28 

37°11 ' 1 2" 

76°03 ' 03" 

CB0109 

Chesapeake 

Bay 

4,12.22,42,72 

371*1 1 '  48" 

76°02 1  or 

CB0110 

Chesapeake 

Bay 

4,17 

37° 13 • on 

76l 01 ' 32" 

CB0202 

Chesapeake 

Bay 

4,12,27 

37°30 '16" 

76  0 1  2 ' 48" 

C30204 

Chesapeal  •' 

Bay 

4,12,32 

37°  39 ' 08" 

/6  OS  8 

CB0206 

Chesapeak 

<jy 

4,22,42 

37°30'15" 

76°Q4 ' 42l 

CB0208 

Chesapeake 

Bay 

4,32,52 

37°30'21 " 

76°02 1 76" 

CB0210 

Chesapeake 

Bay 

3,11,26 

37°30 ’ 28" 

75°59 ' 38" 

CB0301 

Chesapeake 

Bay 

4,22,35 

37°54 ' 02" 

76°  1  1 '42" 

CB0303 

Chesapeake 

Bay 

4,22,62 

.  7°54 ' 1 8" 

76° 1 0 ' 03" 

CB0306 

Chesapeake 

Bay 

4,22,37 

•'7°54’46" 

76°07‘ 16" 

CB0310 

Chesapeake 

Bay 

3,21,57 

3/°56 1 30" 

75°56  1  30" 

CB0401 

Chesapeake 

Bay 

4,22,32 

8l,°23'  14" 

76°22 ' 1 6 " 

CB0404 

Chesapeake 

Bay 

4,42,92 

35°27 1 20" 

76°20 ' 1 3" 

CB0407 

Chesapeake 

Bay 

4,16 

38°23 ' 29" 

76°1 8' 1 5" 

(Cont 

i nued) 

*  Additional  depths 

sampled  during  plan  test  only. 

(Sheet  1  of  6) 


Table  11  (Continued) 


Station 

Body  of  Water 

Sampling  Depths 
(ft  below  msl ) 

Lati tude 

Longitude 

CB0502 

C50505 

Chesapeake  Bay 
Chesapeake  Bay 

4,22,32 

4,52,102 

38°50 ' 32" 
38°50 ' 35” 

76°26'12" 
76°23 ' 50" 

CB0611 

CB0604 

Chesapeake  Bay 
Chesapeake  Bay 

4 ,2S,53 

4,22,37 

39°08' 03" 
39°08' 20" 

76°24'06" 
76°19 ' 23“ 

CB0703 

CB0705 

Chesapeake  Bay 
Chesapeake  Bay 

4,12,29 

3,23,43 

39° 1 9 '57" 
39° 1 9 1 22" 

76°1 2 1 02" 
76° 1 1 ' 24" 

CC0001 

CC0002 

Cape  Henry  Channel 
tape  Henry  Channel 

4,23,52 

4,24,47 

36°59 ' 05" 
36°59 1 48" 

76°00 'll” 
76°00' 54" 

CG0101 

Choptank  River 

4,22,62 

3  8°  3  8 1 10" 

76°09 ' 36' 

COO  101 
C00201 

Chickahominy  River 
Chickahominy  River 

4,20 

4,10 

36°1 6 1 1 3" 
36° 1 8 ’ 42" 

76°52 ' 47" 
76°52 ' 42 " 

EE0101 

EE0201 

EE0301 

E.  Br.  Elizabeth 
E.  Br.  Elizabeth 
E.  Br.  Elizabeth 

Ri  ver 
River 
River 

4,16,28 

4,14  ,24 

4,14,24 

36°50 ' 27" 
36°50 ' 24" 
36°50' 13" 

76°16 ' 56" 
76° 1 6 1  00" 
76°14 1 48" 

EH0101 

EH0102 

EH0103 

Elizabeth  River 

El  izabeth  River 
Elizabeth  River 

4,12 

4,14,24,48,58* 

4,14 

36°55 ’ 30" 
36°55'30" 
36°55 1 30" 

76°20 ' 08" 
76°20  *21" 
76°20 ' 33" 

EH0201 

EH0202 

EH0203 

Elizabeth  River 
Elizabeth  River 
Elizabeth  River 

4 

4,14,24,48,58* 

4,24,43 

36°53'09" 
36°  53 ' 15" 
36°53 1 14" 

76° 1 9 ' 53" 
76°20'04" 
76°20 1 16" 

EH0301 

EH0302 

Elizabeth  River 
Elizabeth  River 

4,14,24.48 

4,14,24 

36°52 ’ 25" 
36°52 ' 23" 

7  6 0 1 9 ' 59" 
76°20 1 12" 

EH0401 

Elizabeth  River 

4,14,24,44,48* 

36°51 ' 31 " 

76°1 8 ' 50" 

EH05Q1 

Elizabeth  River 

4,14,24,42,48* 

36°50 1 34" 

76°17’41" 

EH0601 

Elizabeth  River 

4,14,24,42,48* 

36°48' 50" 

76  °1  7 ' 28" 

EH0701 

Elizabeth  River 

4,14,24,36,43* 

36°47 ' 05" 

76° 18 1 13" 

EH0801 

Elizabeth  River 

4,14,24,36,41* 

36°46 '  36" 

76°1 7 1 45” 

EH0901 

Elizabeth  River 

4,14,24,36 

36° 45’ 18" 

7  6° 1 7 ' 44" 

EH1001 

Elizabeth  River 

3,13 

36°43' 53" 

76° 1 6 ’ 49" 

(Conti nued) 
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Table  14  (Continued) 


Sampling  Depths 


Station 

Body  of  Water 

ft  bel  ow  msl 

Latitude 

Longi tude 

JG0101 

James 

River 

3,13 

36°58'54" 

7  6° 1 7 ‘ 36" 

JG0102 

James 

Ri  ver 

3,23,43 

36°59'36" 

76°17 1 54" 

JG0103 

James 

Ri  ver 

3,13,23,43,72 

37°00' 00" 

7  6° 1 8 ’ 12" 

JG0211 

James 

Ri  ver 

4,14 

36°58' 17" 

76°28'20" 

JG0201 

James 

Ri  ver 

3,12 

36°5  3 ' 48" 

76°27' 54" 

JG0202 

James 

Ri  ver 

3,12,22 

36°58 ' 54" 

76°2  7 ‘ 12" 

JG0203 

James 

Ri  ver 

3,23,43 

36  °5  9 '  1 2" 

76°26 ' 54" 

JG0311 

James 

River 

4,14 

37°03 ' 06" 

76°36 ' 16" 

JG0301 

James 

Ri  ver 

3,12,18 

37°02  54" 

76° 35' 54" 

JG0302 

James 

Ri  ver 

4,13,20,30 

37°03‘ 30" 

76°35 1 36" 

JG0321 

James 

Ri  ver 

4 

37°04 1  11" 

76°36 1 15" 

JG0401 

James 

Ri  ver 

2,19 

37° VI  '  64" 

76°40' 13" 

JG0402 

James 

River 

4,13,28 

3  7° 1 2 ’ 30" 

76°39 ' 12" 

JG0501 

James 

Ri  ver 

3,7,2. 

37°1 2 ' 54" 

76°  48 ’ 0°" 

J  GO  502 

James 

Ri  ver 

3,20,39 

3  7° 1 3 ' 06" 

76°47 ' 30" 

3 GOG 01 

James 

Ri  ver 

3,13,2° 

37° 1 4 ' 06" 

76°56 ' 4P" 

JG0701 

James 

Ri  ver 

4,14,29 

37°17'06" 

77°02‘ 36" 

JG0801 

James 

River 

5,15,25 

3  7° 1 8 ' 06" 

77°08'49" 

JG0901 

James 

Ri  ver 

4,16,29 

37°  3  5 ' 06" 

77°  1 6  1 2 

JG1001 

James 

Ri  ver 

4,13,24 

37°22 ' 48" 

77°20' 18" 

jrioioi 

James 

Ri  ver 

4,12,20 

36°  5  8 ' 06" 

76° 1 9 ' 4°" 

JN  0 102 

James 

Ri  ver 

4,14,24,44,57* 

36°58 ' 1 7" 

76°20' 03" 

JN0103 

James 

River 

4,14,24,34,58 

36  0  58 ' 30" 

76°  20 ' 16" 

J  ,‘10104 

James 

Ri  ver 

4,14,24,44 

36°58 ' 52" 

76°20 ' 41 " 

JN0105 

James 

River 

4,12 

36°59 ' 33" 

76°21 '25" 

JM0201 

James 

River 

4,12 

36°55 ' 09" 

76°25 ' 56" 

JN0202 

James 

Ri  ver 

4,14 

36°55 1 23" 

76°25 1 46" 

JN0203 

James 

Ri  ver 

4,14,23 

36°56 1 24” 

76  25'22" 

JN0204 

James 

River 

4,14,24,50,58* 

36°57;04" 

76°25 1 07" 

JN0205 

James 

River 

4,14,22 

36°57  28" 

76°24 1 41" 

JN0301 

James 

Ri  ver 

4 

37°00  33" 

76°32 ' 44" 

JN03O2 

James 

Ri  ver 

4,18 

37°01  17" 

76°31 ‘ 55" 

JN03O3 

James 

Ri  ver 

4,18,32 

37°01  43" 

76°31 ' 26" 

JN0304 

James 

Ri  ver 

4 

37°02  33" 

76° 30 ' 33" 
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Table  14  (Continued) 


Station 

Body  of  Hater 

Sampling  Depths 
ft  below  msl 

Lati tude 

Longi tude 

JN0401 

James  River 

4,14,24 

37°na 

21" 

76°  39’  53" 

JN0501 

James  River 

4,14 

37°07 

29" 

76°  39'  05" 

JN0502 

James  River 

4,17,31 

37  07 

33" 

76  38'  32" 

JN0601 

James  River 

4,10 

37o10 

30" 

76°  42'  32" 

JN0602 

James  River 

4,10 

37o10 

50" 

76°  43'  10" 

JN0603 

James  River 

4,14,28 

37  11 

15" 

76  43' 53' 

JN0701 

James  River 

4 

37°09 

30" 

76°  44'  02" 

JN0801 

James  River 

4,14,25 

37q1  2 

55" 

Jfo \2' 4!" 

JN0802 

James  River 

4,12 

17°13 

28" 

76052  42" 

JN0803 

James  River 

4,14,24 

37  14 

36" 

76  52' 43" 

LF0101 

Lafayette  River 

4 

0 

36054 

04|; 

76°  19'  11" 

LF0201 

Lafayette  River 

4,16 

36a54 

13 

76nl  7'  40" 

LF0301 

Lafayette  River 

4,10 

36  53 

16 

76  16'  42" 

LH0001 

Lynhaven  Bay 

4.20 

0 

36  55 

09" 

76°05'  18" 

LS0001 

Horseshoe  Shoal 

4,18 

0 

36057, 

26" 

76°  12|  37" 

LS0002 

Horseshoe  Shoal 

4,14,23 

36n58 

53 

7601 2  08' 

LS0003 

Horseshoe  Shoal 

4,14,24 

37o00 

20 

76011'  44" 

LS0004 

Horseshoe  Shoal 

4,18 

37  01 

30 

76  11'  21" 

MB0102 

Mob.iack  Bay 

3,20 

37°  18 

11 

42 

76°  20 '  48" 

MF0001 

Fort  McHenry  Channel 

4,22,52 

39°  1  3' 

n 

48 

76°  32’  32" 

NG0101 

Nanticoke  Rive, 

4,14,24 

38°T^r* 

ti 

50 

75° 55’  39" 

NNOOOl 

Newport  News  Channel 

4,14,32,50,58* 

0 

36  57 

II 

35 

76  21' 16" 

NSOlOl 

Nansemond  River 

4,17 

0  1 

36054, 

12" 

76q27'  32" 

NS0102 

Nansemond  River 

4 

36  54 

23 

76  28’ 16" 

NS0201 

Nansemond  River 

4,12,20 

36°53 ' 

ll 

18 

76°29'  17" 

NS0301 

Nansemond  River 

4,12,19 

36°52 ' 

11 

03 

76°  30 '  3/' 

NS0401 

Nansemond  River 

4,10 

0  1 

36  50 

31 " 

76°  32' 00" 

NS0501 

Nansemond  River 

4 

O  1 

36  49 

ll" 

0  1  „ 

76  32  22 

PGOlOl 

Patuxent  River 

3,22,39 

0  1 

38  18 

43" 

76°25‘  17" 

(Continued) 
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Table  14  (Continued) 


Sampling  Depths 


Station 

Body  of  Water 

( ft  below  msl ) 

Lati tude 

longitude 

P10101 

Piankatank  River 

4 ,12,20 

0  )  n 

37  31  54 

76  °1 8 

’24" 

POO  102 

Potomac  River 

3,22,38 

0 

37058 

'30" 

76°19 

'42" 

POO  103 

Potomac  River 

3,22,40 

37059 

'30" 

76  °1 9 

'54" 

POO  104 

Potomac  River 

3,31,49 

38  00 

52" 

76  °20 

’00" 

P00202 

Potomac  River 

3,22,60 

0 

38  04 

04“ 

76°27 

'll" 

POO  30 1 

Potomac  River 

3,32,57 

0 

38  09 

28" 

76°34 

'39" 

P0040? 

Po t oma c  River 

4,22,41 

0 

38  12 

08" 

76°45 

42 " 

P 00601 

Potomac  River 

3,22,62 

0 

38  21 

20" 

76°59 

'12" 

P00802 

Potomac  River 

5,19 

0 

38  21 

33" 

77°1  4 

59" 

PQ0101 

PoguOSSOii  R i  ve r 

7.10 

0 

37  10 

12" 

76°23 

on" 

PR0103 

Patapsce  Rivet 

1  1  I'D  A<.. 

0 

39  10 

45" 

76°26 

33" 

PR0202 

Pafapsco  River 

1 6,33,50 

0 

39  11 

58" 

o 

ro 

O 

CD 

10" 

RG0101 

Rappahannock  River 

3,13,26 

0 

37034 

54" 

76°017 

36" 

RG0102 

Rappahannock  River 

3,30,33 

37  36 

06" 

76  0  7  7 

06” 

RG0201 

Rappahannock  River 

3,13,26 

o 

37036 

25" 

76°23 

41" 

RG 0202 

Rappahannock  River 

3 , 14,25,52 

37  36 

55" 

76°23 

54" 

RG0301 

Rappahannock  River 

3,13,26,59 

0 

37037, 

54" 

76°28 

36" 

RG0302 

Rappahannock  River 

3,13,20 

37  30 

42" 

76°26 

12" 

RG0501 

Rappahannock  River 

3,13,25 

o 

37  45 

54" 

76°3  7 ' 

12" 

RG0601 

Rappahannock  River 

2,19 

o  , 
37  48 

3b" 

76°42 

36" 

RG0701 

Rappahannock  River 

2,19 

o  , 
37  52 

30  " 

76°46 

18" 

RS 0001 

Rappahannock  Shoal  Ch. 

4,24,52 

o  , 

37034, 

51 " 

76^03 

32" 

RS0002 

Rappahannock  Shoal  Ch. 

4,24,52 

37037, 

00  " 

76^05' 
76  081 

48" 

RS0003 

Rappahannock  Shoal  Ch. 

4,23,52 

37  39 

12" 

02" 

TS0001 

Thimble  Shoal  Channel 

4,14,24,49 

36°56 ’ 

49" 

76°00 

11" 

TS0002 

Thimble  Shoal  Channel 

4,14,24,49,58* 

36°57' 

40" 

76^04 

*6°06 

45" 

TS0003 

Thimble  Shoal  Channel 

4,14,24,49,54* 

36°58 

21" 

’39" 

TS0004 

Thimble  Shoal  Channel 

4,14,24,49,58* 

36^59’ 

37°00 

47' 

/6°1  1 

53" 

TS0005 

Thimble  Shoal  Channel  4,14,24,49,58* 

( Conti nued) 

24"  76  14 1 39" 
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Table  14  (Concluded) 


Sampling  Depths 


Station 

Body  of  Water 

ft  below  msl 

Lati tude 

Longitude 

WB0101 

W.  Br.  Elizabeth  River 

4,14 

36°51 ' 34 " 

76 °20'07" 

WB0201 

W.  Br.  Elizabeth  River 

4,16 

36°51 1 14" 

76°21'20" 

WOO 101 

Willoughby  Bay 

4,10 

36°57‘ 27" 

76  0 1 6 ' 57" 

Y  GO 101 
YG0102 

York  River 

York  River 

4,24,34 

5,25,54 

37°14'12" 
37  1 4 1 30" 

76°25 ' 00" 
76°25 ' 00" 

YG0201 

York  River 

4,14,44,59 

37°  1 4 1 1 8" 

76°30'12" 

YG0301 

YG0302 

York  River 

York  River 

3,13 

3,11,27 

37°19‘00" 
37  19‘18" 

76°36 ' 24" 
76°35 ' 54" 

YG0401 

YG0402 

York  River 

York  River 

3,11 

12,18,25 

37°22‘ 12" 
37  22*24" 

76°39 ' 00" 
76°38 ' 36" 

YG0501 

York  River 

4,14,26 

37°24 1 48" 

76°41 '06" 

YG0601 

York  River 

4,14,24 

37°29 ' 1 8" 

76°45'24" 

YS0001 

YS0002 

YS0003 

YS0004 

YS0005 

York  Spit  Channel 

York  Spit  Channel 

York  Spit  Channel 

York  Spit  Channel 

York  Spit  Channel 

4.24.49 

4.22.49 

4,22,49 

4,22,49 

4,24,52 

37°02'33" 
3r06'24" 
37^1 1 ‘ 03" 

3  71’  1 4 1 12" 
37° 1 7' 12" 

76^04 '31" 
76°07' 16" 
76°09 ' 1 8" 
76°08' 08" 
76°06 ' 59" 
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